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Abstract. We review some recent developments in nonperturbative studies 
of quantum field theory (QFT) using the Schwinger-Dyson equations formu- 
lated directly in Minkowski space. We begin with the introduction of essential 
ideas of the integral representation in QFT and a discussion of renormaliza- 
tion in this approach. The technique based on the integral representation of 
Green's functions is exploited to solve Schwinger-Dyson equations in several 
quantum field models, eg. in scalar models and in strong coupling QED^^i in 
the quenched and in the unquenched approximation. The phenomenon of dy- 
namical chiral symmetry breaking in regularized theory is touched. In QCD, 
the analyticity of gluon propagator on the complex momentum square plane 
is exploited to continue some recent lattice data to timelike momentum axis. 
We find non-positive absorptive part contribution in the Landau gauge gluon 
propagator which is in agreement with some other new recent analyzes. 



1 Introduction 

This work is devoted to studies of quantum field models with strong coupling. 
The Schwinger-Dyson equations (SDEs) in momentum representation are solved 
and analyzed in Minkowski space. In this introduction our motivations, broad 
contexts and main results are briefly summarized. 

These days the experimentally accessible particle physics is being described 
in the framework of the (effective) quantum field theory (QFT) by the Stan- 
dard model of the strong and electroweak interactions with the gauge group 
SUc{3) (8) SUl{'2) ® Uy{X)- When the perturbation theory is applicable the ob- 
servables follow from the Lagrangian in relatively straightforward way: the corre- 
sponding S-matrix elements are generated from n-point Green's functions (GFs) 
which are calculated and regularized/renormalized order by order according to 
the standard rules. Order by order perturbation theory can be naturally gener- 
ated from expansion of SDEs: an infinite tower of coupled equations connecting 
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successively higher and higher points GFs. The GFs themselves are not ob- 
servable and they contain some redundant information: off-shell behavior (con- 
nected e.g. with a possible field redefinition) and also some other remnants of 
the particular calculation scheme employed, e.g., gauge fixing parameters and 
renormalization scales. If the program of the perturbation theory is carried out 
consistently up to the certain order, the observables should be independent of 
all these unphysical parameters. 

Needless to say, many phenomena exist outside the reach of the perturba- 
tion theory approach, especially when the strong sector of the Standard model - 
Quantum Chromodynamics (QCD) - or some extensions of the Standard model 
are concerned: bound states, confinement, dynamical mass generation etc. Going 
beyond perturbation theory requires development of elaborated approaches and 
sophisticated tools, especially in the strong coupling regime. Number of such ap- 
proaches has been pursued: lattice theories, Feynman-Schwinger representation, 
nonperturbative treatment of SDEs, bag models .... These days neither of them 
provides such a clear and unambiguous understanding of the physics as we are 
accustomed to in the perturbative regime. The lattice approach clearly stands 
by its own, since it does not rely on any uncontrollable approximations: just the 
brute force discretization of path integrals is employed. Lattice "data" should (in 
principle) provide as good representation of the studied dynamics as real ones 
(with additional merit of allowing to vary the dynamical input), but converged 
unquenched results are in most cases just not available as yet (however for some 
recent progress in QCD see [1] , |2] , [B] ) • 

In this work we deal with the nonperturbative solutions of the SDEs. If 
one could solve the full infinite set of SDEs, the complete solution of the QFT 
will be available, coinciding with the results of converged lattice calculations. In 
reality one has to truncate this set of equations: the main weakness of SDEs phe- 
nomenology is the necessity to employ some Ansatze for higher Greens functions. 
The reliability of these approximations can be estimated by comparison with the 
known perturbation theory result in the regime of a soft coupling constant and 
by comparison with the lattice data (and/or results of alternative approaches 
like (^-derivable effective action, etc.) where available. The recent results (see 
e.g. reviews p, 5, 6, 7 and references therein) provide some encouragement: they 
suggest that SDEs are viable tool for obtaining (eventually) the truly nonpertur- 
bative answers for plethora of fundamental questions: dynamical chiral symmetry 
breaking, confinement of colored objects in QCD, the high temperature super- 
conductivity in condensed matter (modeled |S] by QED2+i)- In the QCD sector 
there is also a decent agreement with some experimental data: for meson decay 
constants like F^^ , evolution of quark masses to their phenomenologically known 
(constituent) values, bound states properties obtained from Bethe-Salpeter or 
Fadeev equations (the part of SDEs). 

The recent solutions of SDEs 0,12] nicely agree with the lattice data 
[21 El ^3 (the explicit comparison is made in |U1 [71 Ej). Thus, encouraging 
connections between results of SDEs studies and fundamental theory (QCD) 
and/or various phenomenological approaches (chiral perturbation theory, con- 
stituent quark model, vector meson dominance . . . ) are emerging. 
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Most of these results were obtained by solution of SDEs in Euclidean space. 
One of the main goal of this work is to develop an alternative approach allowing 
to solve them directly in Minkowski space. In this method (generalized) spectral 
decompositions of the GFs are employed, based on their analytical properties. 
The main merit of this approach is possibility to get solutions both for spacelike 
and timelike momenta. Techniques of solving SDEs directly in Minkowski space 
are much less developed than the corresponding Euclidean ones. Therefore most 
of this paper deals with developing and testing such solutions on some simple 
QFT models. 

For large class of nonasymptotically free theories running couplings - calcu- 
lated from perturbation theory - diverge at some finite but usually large spacelike 
scale. If perturbation theory results are taken literally, stretching them clearly 
beyond the range of applicability, also the QCD coupling develops singularity for 
low momenta squared. In the latter case the lattice as well as SDEs results sug- 
gest that nonperturbative effects alter behavior of running coupling oiqcd{<1^) 
significantly, situation for nonasymptotically free is less clear. This is why we 
pay special attention to behavior of our solutions in the region of high momenta. 

The paper is organized as follows. Section 2 serves as an introduction: some 
well known general results and useful textbooks formulas are collected and basics 
of the formalism employed later are discussed. The integral representation for 
GFs are reviewed and the calculation of dispersion relation is exhibited at one 
and two loop example. Section 3 is devoted to the simple scalar model + 
for pedagogical reason. The momentum space SDE is converted to the Unitary 
Equation for propagator spectral density. The Unitary Equation is shown to be 
real equation where the singularity accompanying usual Minkowski space calcu- 
lation are avoided. The relativistic problem of bounds states as an intrinsically 
nonperturbative phenomena is briefly reviewed in this Section too. 

Section 4 represents study of the SDEs in 3+1 QED. For this purpose we 
improve and reorganize the Unitary Equations already particularly found in the 
papers jl3l I14j in a way that the numerical obstacles with singular integral ker- 
nels are avoided. Due to this the resulting numerical solutions of our Unitary 
Equations become more stable and accurate. 

One should mention at this point a possibility of getting the solutions of 
the SDEs in Minkowski space with the help of analytical continuation of the 
Euclidean ones. Nevertheless in simple cases it can be successful: we have recently 
done this for the gluon form factor and we present our results in the Section 
devoted to the QCD. 

Section 5 offers some further comments on spectral treatment of Schwinger- 
Dyson equations and possible further progress is discussed. 

Section 6 presents the known attempt to solve the SDEs in QCD in Minkowski 
space. Furthermore, S'C/(3) Yang-Mills, i.e the pure gluodynamics is discussed 
there and analytical behavior of gluon propagator is deduced by the analytical 
continuation from Euclidean space to the analyticity domain boundary. To obtain 
numerically stable analytical continuation is not an easy task and attempts to 
perform it are rare. As an input data we have used recently obtained lattice data 

in- 
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2 Integral Representations in Quantum Field Theory 

Solutions of the SDEs in this work are based on integral representations of 
Green's functions. In this Section we briefly review their derivation and then 
present two illustrative examples from the model. 

The development of spectral techniques can be stimulated by following con- 
sideration. Let us assume that any GF can be written as an expansion of skeleton 
contributions, i.e., the diagrams where the propagators corresponding to the in- 
ternal lines are exact. If we further assume the analyticity of these propagators 
( in the sense discussed in the next text), then one can apply the whole well 
established knowledge about the analytical structure of Feynman diagrams. The 
singularity structure of each Feynman graph allows a simple connection between 
momentum dependence in spacelike and timelike regime, i.e., between Euclidean 
and Minkowski space. That the equations of motions which hold between Green's 
functions, i.e. SDEs , can be translated into the corresponding relations between 
spectral functions was already known half century ago jl5l I16j . The renormal- 
ization program has been carried in the terms of spectral function as well. 

This was also the basis of historical development of the dispersion relation 
(for summary see [HJ and of the Mandelstam representation in the late 50's |18j . 
Many of these efforts have been summarized in the standard textbook of that 
time, see e.g. [19'^. This book represents deductive approach to S-matrix theory, 
the subject rather popular in particle physics at that time. Fortunately, being 
mainly based on the analysis of Feynman graphs (and their infinite sum ), the 
methods reviewed in this book ( including also asymptotic high energy behavior, 
Regge poles,..) survive and have found their appropriate use, even after the time 
when QCD was born and the trust to quantum field theory was renewed. 



2.1 The Kallen-Lehmann Representation 

In this Section we remind the derivation of the Umezawa-Kamefuchi-Kallen- 
Lehmann spectral representation (SR) j2()| I2H \22\ l2,Sj for a two-point Green's 
function. 

Let us consider the vacuum expectation value of the product of two real scalar 
fields < 0|(/>(x)(/>(y)|0 >. To arrive at the desired expression for the propagator 

iG{x - y) =< 0\T^{x)ct){y)\0 > (1) 

one inserts the identity operator 

i = |0 >< 0| + y f |Ap >< Apl , (2) 

' A (27r)3 2Vp2+m2' ^ ^' 

between the fields (/)(x) and (j){y) in (|2.H1 : 
< O|0(x)0(y)|O >=Y. f ,,fL . < 0|</'(:r)|Ap X Xp\<P{y)\0 > (3) 

where the sum runs over complete set of A-particles states |Ap > with given 
momentum p. It is assumed that a spontaneous symmetry breaking does not take 
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place, i.e., that < 0|(/>(x)|0 >= 0; if it does a space-time independent constant 
appears on the r.h.s. Making use of the translational invariance, transformation 
properties in respect to the Lorenz boost: 

< 0\(t>{x)\\p >=< 0|e'^-^(A(0)e-^-^-^|Ap > 

=< 0|<A(0)|Ap > e-^npo=E, =< 0|<A(0)|Ao > e-'npo=E, (4) 



and of the integral representation of the step function, one gets for xq > yo 

J ^^^> f ~ '""A 

(5) 

Taking together with an analogous expression for an opposite time ordering 
gives the Feynman propagator in the form 



< 0|^(.)^(,)|0 ^ E / "'J/^J'^J^^ l< 0|^(0)|Ao >| 



G{x — y) = J du) a{u>)D{x — y,Lo) , (6) 
ctH = E5(u;-mi)|<O|0(O)|Ao >|' , (7) 



A 



G{P) = [ du , . , (8) 

J p-^ — uj + ie 

where ^(a;) is a positive spectral density {the Lehmann function) and D stands 
for free-particle propagator with 'mass' -y^iJ,i.e. 

r,f ^ f d^P exp{-ip- {x-y)) 

D{x-y,uj)= ——^ 2 ^) 

J (2vrj^ p"^ — uj + ze 

and D{p,uj) is its Fourier transform. 

To sum it up: the derivation of the spectral representation of the two-point 
Green's function is based on three assumptions: 

• the Poincare invariance of the theory, 

• existence of the physical states |Ap >, 

• the completeness relation 0. 

The physical propagator usually contains the single particle pole and the contin- 
uous part stemming from the interaction. Therefore we typically get: 

G{p) = j duja{u)D{p,uj) =rD{p,n?)+ j duja{uj)D{p,uj) , (10) 

where D{jp,m?) is a free propagator with the physical mass, r is a residuum 
at p^ = m?, and cr is a nonsingular continuous part, nonzero above certain 
kinematical threshold. For the free scalar field r = 1 and £7 = 0. 




Figure 1. The position of singularities on the complex momentum plane with the sketched 
integration contour C used in the main text. The singularities can appear only for Rek^ > 
(crosses) . 



2.2 Analyticity 

Suppose that the complex function f{z) is analytical in the whole complex plane 
of z = fc^ except for the positive real axis z > 0. Then we can choose a closed 
loop C in the complex k'^ plane as in Fig. ^ such that the complex function f{z) 
is analytical inside and on the loop C . For a reference point k'^ inside C, the 
Cauchy integral formula tells us that 



1 f"" ^ J{z + ie)-f{z-ie) ^ 1 /■ f{z) 
2vri Jq z-k"^ 2m J\z\=r z - k"^ ' 



where we have separated the integral into two pieces: one is the contribution from 
the paths above and below the positive real axis; the other is the contribution 
to the integral over the circle of radius R. 

First, making use of the fact that f{k'^) is real for k'^ < Smin for some nonneg- 
ative Smin (by assumption, at least for k'^ < 0, i.e., for spacelike A;), it follows from 
the Schwartz reflection principle that for Re z > Smin 

it holds /(z-ie) = f*{z+ie) 
and hence f{z + ie) — f{z — ie) = f{z + ie) — f*{z + ie) = 2i Im f{z + ie). Hence 
the first term in r.h.s. of H12|) picks up the discontinuity (the imaginary part) of 
the function f{z) along the positive real axis. 

Next, the circle \z\ = R is parameterized as z = Re^^: 



R 



, Im f(z + ie) 
z-k^ + 



2^ de Re'^f{Re^^) 
2^ 



id 



k^ 



(13) 
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d9 Re^'fjRe^'^) f'^ dO R fjRe^ „ 



R-\k'^\ ■ 
(14) 

Now we wish to take the radius R of the circle to infinity. Thus, if |/(^)| — > as 
1^1 — > oo, then the second term has no contribution as ii ^ oo and we arrive at 
the desired Dispersion Relation (DR): 

2, 1 r... 



fin = - dz (15) 



This can also be written as 



oo 



/(^')=/ ds , pis):=-lmfis), (16) 

Jo s - k^ + ie IT 

then the SR for propagator (|lflj) exactly coincides with the DR (|16j) . 

To solve the SDEs in Minkowski space we will need also integral representa- 
tion for the inverse of propagator. In the renormalizable theories such spectral 
functions often do not vanish at infinity. This can be already viewed from the 
identity 

lmf{z)-' = -\f{z)\-'lmf{z), (17) 

where we identified the function / with the propagator function G. Then one 
makes the momentum space subtraction(s) and arrives at modified dispersion 
relations. For instance, let us define a new function 

fB.{C;z) = f{z)-fic), (18) 

and assume that fniz) 0. This gives the once-subtracted DR: 

If the other divergence occurs, one can go on and subtract the next term in the 
Taylor expansion of f{z). 



2.3 Dispersion relations in scalar theory 

In this Section we derive the DRs first for the one-loop elastic two particle scatter- 
ing amplitude in the scalar (p^ theory without spontaneous symmetry breaking, 
i.e., in the QFT given by Lagrangian: 

C = ^5^^a^<?- V(^), (20) 

v(^) = !li<p^ + go<p\ 

where the global minimum condition = 0) = entails the zero value of the 
field vacuum expectation value. 
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Next, as a less trivial example, we offer a detailed original proof of DR for 
a two-loop self-energy function in the (off-shell) momentum subtraction renor- 
malization scheme. The calculation is performed for the noninteger dimension 
D = 4: — e, the limit e — > is taken after the appropriate subtractions. However 
the results on imaginary part of the sunset diagram was already obtained several 
years ago by the Cutkosky rules approach 24 (see also |25i for general case in 
n-dimensions) the derivation presented here follows an ideas of renormalization 
realized directly in momentum space. 

2.4 Scattering amplitude at one skeleton loop 

The scattering matrix in theory is composed from the four-point Green's 
function. The one-loop skeleton contribution to the irreducible GF reads 

.(4),,2^_ V V _ .. (4!go)^ f d^k 

{2ttY 




n 'il')= I X =-^'-^ I 7KZr^G{k + l)G{k), (21) 



where go is a bare coupling and G represents fully dressed propagator (jlUj) . 
Substituting the Lehmann representation for G gives 

= ^^j^jdaiaiai) Ida2a{o^2)lo{l'), 



m') = I ^Dik + l,a,)D{k,a2). (22) 



Thus, the derivation of DR reduces to the familiar "perturbative" spectral rep- 
resentation of the scalar integral /o(^^)- The unrenormalized vertex Fq requires 
just one subtraction due to the presence of logarithmic UV divergence 

r^^Hl^Cg) = ^^^^^ j daia{ai) j da2a{a2)lR{l^), 



QO 



= / ^(l^lJf^-h^ ,23) 

J {uj - Cg)[l^ - u + le) 

where the function X is defined as 

X(a, b- a;) = ^'^'^"'^'^^ 0{u - + Vbf) , 

X{x, y, z) = x^ + + - 2xy - 2yz - 2xz . (24) 

and (^g is the (square of) renormalization scale. In order to define the renormal- 
ized coupling g{Cg) we choose the function r^^\PXg) to be renormalized loop 
contribution of quartic vertex and include the infinite constant -ro(Cg) into the 
counter-term part of go. The full scattering amplitude M(s, t, u) in the considered 
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approximation is the sum of the infinite chain of the irreducible contributions 
defined above: 

Mis,t,u) = Alg{C,)+ Y^^^^t^y (25) 

x=s,t,u ^ ' ^5/ 

where one common renormafization scale Q was used to renormalize GF. F in 
each channel characterized by the Mandelstam variable s, t, u (labeling incoming 
(outgoing) momenta pi^2 (P3,4)) they are s = (ipi +p2Y-,t = {pi —p^Y,u = 
{pi — Pif')- The first term on the l.h.s. of Eq. (|25j) represents the tree value of 
the amplitudes M. 

Several important remarks can be made about the obtained result. 

• Integrating the single particle contribution we explicitly see that the vertex 
becomes complex for the timelike x > 4m^ and its absorptive part reads: 



167r 



/ f°° 
r^Wl ©(w - 4m^) + 2r / da a{a)X{a,m'^ , 



+ / da2cr{a2) / dai a{ai)X{ai,a2,u>) 

which is explicitly independent on the renormalization scale Qg. 

• perturbation theory is trivially reproduced. Assuming r ~ 1 the one-loop 
result in the perturbation theory is given by the first term in (|26() . 

• The absorptive part does not vanish for large momenta, it is a constant. For 
fixed momentum the naive limit C,g ^ oo cannot be performed . Taking this 
explicitly then the original log. divergence reappear again. It is impossible 
to define infinite momentum frame renormalization scheme (the so called 
Weinberg scheme j26j). 

• The summation of the irreducible contributions leads to the Landau ghosts 
in the amplitudes and the renormalized coupling can not be reasonably 
defined without some regulators. This pathology is manifest already at the 
one-loop level. This is because the theory is the (probably the best 
known) example of the nonasymptotically free theory. 

2.5 DR for sunset diagram in momentum subtraction renormalization scheme 

In the relation (|2(ij) knowledge of the Lehmann function a is required. It can be 
taken from the perturbation theory or some improved approximation of the full 
result can be used. To this end we continue with derivation of a similar relation 
for the two-point Green function. The renormalized inverse propagator is given 
by 

G~\p')=p'-m\C)-nR{C;p^), (27) 
where is the renormalized two-point proper self-energy function: 



nn{C;p^) = n{p')-n{0 



dn{p 



2\ 



dp^ 



ip' - C) . (28) 



10 



Implications of analyticity to. 



since the quadratically divergent self-energy is still logarithmically divergent even 
after the "mass subtraction" . This 'infinity' proportional to is to be absorbed 
by the field-strength redefinition. The two loop proper selfenergy function to be 
renormalized is: 



Here, it is taken in the skeleton sunset approximation: the corrections to the 
vertex are neglected and the full quartic is replaced by its tree value -Tj^ee = 'ilg- 
The direct consequence of the subtraction procedure (|28|) is the DR formula 
for the renormalized self-energy 

nMC^p')= f^^Jl^^Jp:^. (30) 

Our task here is to find the weight function p{uj) for the skeleton sunset approx- 
imation considered. As in the previous one-loop case it is convenient to drop all 
prefactors and thus to perform effectively the subtraction (|28|) at the level of 
"bare propagators" : 

no{p^) = - ^^'^^ j daia{ai) j da2a{a2) j da^aias) J{p^) , 

/d'^k f d'^k 

(31) 

Our goal is to cast the function J{p'^) into the integral over some spectral vari- 
ables, i.e., to manipulate integrals so that integrations over momenta can be 
taken. While the one loop case H23|) proceed standard textbook derivation (note, 
the details of the derivation are provided in the Section 3, follow the relations 
H54|l to ()58() ) the derivation the two- loop derivation is less trivial. We do it for 
this particular example in detail in the Appendix A, here we present the result 
firstly published in [23|,[1S]- The absorptive part of the integral J{uj) reads 

(v^-v^) / 3 \ 

/ X 1 f J \/A(s,u;,a3)v'A(a2,s,ai) / 1/2.3 \ 

(32) 

which leads to the full expression for the spectral function in DR (|Hn|l for sunset 
skeleton graph: 

Im77K(p2) (4!^) 



vr 



^^I^YlJdaMai) 




VA(.,p^a3)A(«2,. ^^ j ^-2 _ ^^^y j _ (33) 
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In massive t/)^ the spectral function is nonzero when > 9m^, otherwise it is 
zero. 

2.6 Perturbation Theory Integral Representation in renormalizable theories 

The Perturbation Theory Integral Representation (PTIR) |221 is a natural ex- 
tension of the Lehmann representation for the propagator to n— point Green's 
functions. The PTIR allows to express the amplitudes as integrals over weight 
functions and a momenta-dependent denominator with known singularity struc- 
ture. Since the Feynman parametric integral always exists for any Feynman di- 
agram defined by the perturbation theory, one can always define an integral 
representation such that the number of independent integration parameters is 
equal to that of invariant squares of external momenta. The general expression 
for n-point amplitudes, which is considered to be a sum of all Feynman diagrams 
with n external legs, has been derived by Nakanishi [JJj for scalar theories with 
nonderivative couplings. The PTIR is unique for unrenormalized amplitudes. The 
general proof of the Uniqueness Theorem can be also found in the Nakanishi's 
book The renormalized 2,3,4-point vertices can differ only by the presence of 
subtraction polynomials as follows from renormalization procedure. The coeffi- 
cients of the subtraction polynomials then depend on the renormalization scheme 
employed. 

For theories with interactions containing field derivatives or theories with 
particles of nonzero spin similar results to all order of perturbation theory have 
never been proved, since tensor and spinor structure of corresponding expressions 
becomes rather complicated with increasing number of external legs of Feynman 
diagrams. We do not feel too discouraged by this: the (generalized) PTIR can 
always be deduced for a given Feynman diagram. Indeed, if we assume the an- 
alyticity of the lowest point Green's function -the propagators- then it is sure 
that the PTIR is derivable for skeleton diagram under consideration. 

In what follows we briefly review examples of PTIR for scalar amplitudes. We 
start with PTIR for the two-point correlation function 11 (p^), already discussed 
in detail above. Recall that the expression for unrenormalized self-energy reads 



where the finite coefficients a, b are characteristic of a given but arbitrary renor- 
malization scheme. Their size indicates how a given scheme employed differs 
from the momentum subtraction renormalization scheme where a,b = 0, since 
the coefficients 




(34) 



The renormalized UrIp^) turns out to be 




(35) 




(36) 



are fully absorbed into the Lagrangian counter-terms in this scheme. 
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For the three- leg vertex function r{pi,p2) analysis of contributing Feynman 
diagrams leads to the PTIR derived by Nakanishi: 



3 3 



P3{a,z) 



a - {zipj + Z2pI + Z3pI) - ie 



(37) 



where the momenta are conserved Pi + P2 + P3 = 0, the invariant squares are 
independent and the single weight function p3{a,z) is sufficient to describe the 
sum of all relevant Feynman diagrams. 

Further analysis of this vertex function gives various constrains on the inte- 
gration regions a,z, dependent on masses of the particles. For instance, for the 
^^(p interaction and m (//) being the mass of particle ^ (i;^), Eq. H37() reduces to 
the following two-variable spectral representation: 



r{p, P)= da dz 



p{a,z) 



a 



-1 



(p + z^y 



(38) 



in the region < < 4m^, P = pi + P2i P = Pi — P2- The spectral function 
p{a, z) is a positive regular function of the spectral variables for the fixed value 
of P^. Furthermore, if the internal lines with momenta pi and p2 correspond 
to the particle <P and the third external line corresponds to the particle the 
following support of p{a, z) can be derived: 



p(a, z) = 0, unless a> \ m + p — {1 — \z\) 



(39) 



This two variable representation was first proposed in an axiomatic field theory 
[2Hj,[2ni and successfully used to solve the Bethe-Salpeter equation in various 
approximations |5n],|3Tj. The paper solves both the bound state equation 
and the scattering problem and applies this two variable representation to QED 
even in the nonperturbative regime. 

The scattering matrix M(p,q;P) describes the process ^1^2 ^3^4, where 
p = q\ — q2 and q = q^ — q^ are the initial and final relative momenta, respectively, 
and P is the total four-momentum. It is given by the infinite series of Feynman 
diagrams. The full renormalized scattering matrix can formally be written as 



M{p,q;P) 



+ 



n I 



le 



Pus{i,i) 



(40) 
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where is the square of the four-momentum carried by 'Pi and s, t, and u are 
the Mandelstam variables. Since only six of these variables are independent (due 
to the relation ^i+^l + gl + ^l = s + t + u), this is also the number of independent 
^—parameters. Hence, one only has to introduce the "mass" parameter 7 and six 
dimensionless Feynman parameters with one constraint. The symbol 17 denotes 
the integration region of ^ such that i7 = {^| < .^j < 1, ^ = l(i = 1, ...,6)}. 
The function pch gives the weight of the spectrum arising from the different 
channels which can be denoted by ch = {st}, {iu}, {us}. Since any Feynman 
diagram for the scattering matrix can be written in this form, it must be also true 
for their sum. From the corresponding SDEs of considered theory it simply follows 
that PTIR for M must contain the weights of lower point Green's functions. For 
instance, the simplest contribution from f-channel is given by a free propagator 

M{p,q;P) = — f . , (41) 

Kxch -{p-qY-ie 

and the corresponding weight is simply a product of six delta functions. 

The first topologically less trivial diagrams are various boxes and crossed- 
boxes. Collecting various couplings to the constant prefactor C then the two- 
variable expression for crossed box diagram was derived some time ago by Man- 
delstam [H 



C 



00 poo 



dt' du' e{K{u',t')) 



74^2 Jim2 t' -t-ieu' -u-ie ^ k{u' ,t') 
K{u',t') = u't' [{t' - Afi^){u' - 4m2) - Afi^] , (42) 

where fi is mass of the exchanged particle (f) (its propagators cross inside the 
diagram), while m is a mass of the scattered particle (emitting and absorbing 
(j)). Here, the function 

a{u,t) = C , , , , (43) 

is called the double spectral function. It can be easily converted to the PTIR 
form: 

00 1 

Icross = [da [dz . , (44) 

J J a — zu — [1 — z)t 


with the help of the Feynman identity. 



3 SDEs and Unitary Equations in the scalar toy model 

This Section should serve as a pedagogical introduction to the more general 
strategy of Unitary Equations. As we will see, using the spectral decomposition 
of GFs in SDEs allows to derive a real nonsingular integral equation for a real 
spectral function a. For this purpose we consider a simple quantum model de- 
scribing the self-interaction of one real spinless field (p. The Lagrangian density 
for our toy model reads 

jC = ^d^Mx)d^'Mx) - lrnl4{x) - go4{x) + Xo4{x) + ... , (45) 
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where the subscript indicates the unrenormalized quantities and the dots mean 
the neglected interactions with neglected field content. 

The main philosophy of the derivation does not depend on the details of 
a given model. The method can be directly applied to all theories where the 
dispersion relation approach is accessible. We also explain and review the spectral 
technique used in the calculations of relativistic bound states in the subsequent 
Section. For this purpose we follow the main ideas of the paper [H^ where the 
(massive and gauged) Wick-Cutkosky model and its response on electromagnetic 
interaction was considered. Furthermore, we not only exhibit the appropriate 
derivation of Unitary Equations but we also present the numerical solution and 
we offer the comparison with the results obtained via usual Euclidean formalism. 

Using the diagrammatic representation for the selfenergy the SDE for scalar 
propagator G can be depicted as: 

Go'(p')=/-"^o-^o(/), (46) 



no{q')= ^ + 



(47) 

where the full irreducible vertex functions F^"^ i^t^) ^-^e represented by blobs 
(diamond) and they satisfy their own SDEs (for the derivation of SDEs form 
the effective action see some standard textbook, e.g. [SSI, but also [Sll)- All the 
internal lines stand for G ,i.e. for fully dressed propagators. For the explanatory 
simplicity we neglect the two loops diagram, how to take into account the last two 
diagrams will be explained at the end of this Section. In our simple approximation 
the selfenergy to be considered is given by the first two diagrams above. Their 
contribution is: 

^,{v')=^\ I ^Go(g)+i3<7o J ^,rl,'\p - q,q)Go{p - q)Go{q) . (48) 

Furthermore, we approximate the trilinear and quartic vertices by their tree 
value: -T'^I = 3!^ and F^^'^ = 4!A (note, contrary to the usual convention used 
in QED, we very naturally include the coupling constants into the definition of 
proper GFs). In such approximation the quartic interaction \<P'^ as well as the 
field do not need renormalization since we have neglected higher diagrams. Thus 
the divergent integral defining TTq requires only the mass renormalization: 

nniCp^) = noip^) - noiO , (49) 

where the divergence (quadratic here) is absorbed to the counterterm part of the 
Lagrangian: 

ml = m2(C) + Sm^ , 6m^ = 77o(C). (50) 
The renormalized equation for the propagator then reads 

G-\p^) = p^-miCf-nRiC;p^), 
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nRic-^p") = n{p^)-n{C), (5i) 
n{p') = i2,g j ^^r^Mp - q)G{q) , (52) 

where we have omitted the seagull diagram contribution in expression for 11 {p^) 
since it exactly vanishes after the subtraction and where we simply take Go(p^) = 
Gip^) since the field strength renormalization constant Z = 1 here (for more 
complete discussion of the multiplicative renormalization in scalar models see 
the paper [35]). 

The key point of our spectral technique is the derivation of the Unitary Equa- 
tions. The Unitary Equations relate the imaginary part of the propagator -the 
spectral function a- with the imaginary and the real part of proper GF (this 
is the reason why we call them the Unitary Equations). Dealing with explic- 
itly massive model here, the generic spectral decomposition of the renormalized 
propagator ((TU)) is used, where the spectral function of the stable particle is 

a{a) = r5{m?' — a) + a{a) . (53) 

When propagator SR is used in the equation for 77 Eq. (|52|) it leads to the 
DR for the renormalized selfenergy nii{C,;p'^). Its form is well known from the 
perturbation theory treatment and is given by the DR already discussed in the 
Section 3. However here, for purpose of convenience of the reader and in order 
to make our explanation self-contained, we also give the appropriate derivation 
here. 

Substituting SR for propagators into the selfenergy U reads 

i7(p2) = l8c?2 j dadl3a{a)a{l3)I{p^) , (54) 

^ J (27r)4 {{p + /)2 - a + ie){P - p + ie) ' 

Clearly nji{(;p'^) is given by the integral over the analogically defined Ir{C'iP'^) = 
I{p^) — /(C)- Using the Feynman parameterization Ir can be written as 



1 dH_ 

(27r)4 1^ [/2 _^ p2^(^i -x)- a(l -x)- I3x + ie] 



UC;p') = I dxjj^^{ . , ,^ ; . _ . . - (p^ - C) 



1 

dH f —2idxdz(p^ — C)x{l — x) 



(27r)4 J [p + (p2 _ Q^i^i _ j_ ^^(1 -x)- a(l -x)- I3x + ie 



3 







^.^i^f^, (55) 

Z — II + l£ 

where we have labeled fl = C, -\ — C+<^/x+^P/{i x) ^ jyjg^j^jj^g substitution z — > u; 
such that uj = Q, \i yields 

I rl f-oo p^ — C 

UC;p') = j^l '^^A,^ (,2 + (56) 

X 1 — x 
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Changing the ordering of the integrations (note, the variables a, (5 are positive) 
and putting / into 77 we obtain the desired DR: 



duj — 7^ — 







■ u) + ie 



(57) 



where n = 1 here (while n = 2 when the field renormalization enters the calcula- 
tion and the second subtraction is required, here we use n for general purpose) 
and where the absorptive part up is given as 



p{uj) 



aaapa{a)a{p) &{lu 



ijj 



(^/^+a/^)2) 



4m ^ 



0(a; - \n?) + . . . , 



(58) 



where A is the Khallen triangle function defined in ()24() and the dots means 
the terms with one or two it's. (Thorough this paper we always use the letter 
p to denote the weight function that appears in selfenergy in order to carefully 
distinguish from the spectral function a , i.e. from ImG(p^)/7r). Cg is used to 
label the coupling strength defined as Cg — 



'9 — JA^- 

The derivation of Unitary Equations is achieved by using the functional iden- 



tity 



1 



P. 



1 



iit5{x' — x) . 



(59) 



X — X + le X — X 
where P- denotes the principal value integration. GG~^ = 1 can be written as: 

n 



p2 _ ^2 



+ 



da 



a — le 



— m^((") 



duj- 



p{uj) 



(jj + ie 



1 



(60) 

we arrive at the equations between the spectral functions a and p and their 
principal values integral. We eliminate the unpleasant (since unknown) principal 
value integral over the function a from the equation for the real part 



P. / did 



a{u;) l-7r^pip^)a{p^) 



UJ 



h{p^) 



(61) 



h{p') =ReG~^(p 



p 



n?{C) -P. I duj 



p{u) 



9 

p'^ — UJ 



The equation for imaginary part of Eq. (|60|) can be written as 



a{p')b{p') - p{p') 



+ P. dio 



p^ — UJ 



(62) 



Substituting (|HT|l into the Eq. H62() we arrive at the desired Unitary Equation: 
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as follows from the rel. (|58)1 . the continuous spectral function a starts to be 
smoothly nonzero from the perturbative threshold Am?. Note that the function 
b still contain the principal value integral this can evaluated numerically |SS| 
igeneral, but in simple one loop case here it is known analytically. The form of 
the Eq. (|63|) is based only on the analytic property of the propagator, it does not 
depend on a details of interaction. All the dynamical information is contained in 
the expression for p which follows from the SDEs solution here, or from the usual 
Feynman diagram in the case of perturbation theory treatment. At the end we 
should mention that the derivation of Unitary Equations has been performed for 
momentum subtraction scheme only for simplicity. It should be stressed here, 
that the extension to any other renormalization scheme is very straightforward. 
Remind the changes in DR which arise due to finite change in the field renormal- 
ization constant Z. It would be reflected by the presence of nontrivial subtracting 
polynom before the DR This can be easily incorporated into the Unitary 
Equations. 

If the pole is assumed in the propagator then the knowledge of the pole mass 
m is required in the presented spectral treatment. From its deflnition G~^(m^) = 
we can obtain the following equation 

J — I 

The propagator residuum can be obtained by the on-shell differentiation: 

2 2 -I 

p — m 1 
T — lim = 

J iuj—rrf')^ 



(64) 



3.1 Comparison with the Euchdean solution 

It is instructive to compare our Minkowski-space result with the one obtained in 
Euclidean approach. For this purpose we consider the propagator SDE truncated 
and renormalized as before. 

After the Wick rotation and angular (3, (p integrations of the four-dimensional 
sphere 



1 

2^ 



desin^iO) / d/3sin(/?) 



27r 



(66) 



the renormalized SDE for propagator becomes 



7T(x) = 771^(0) 



2Ca 



IT 



dy / d0 



sin^(0)y 

y + n[y) 



1 



z + n{z) y + n{y) 



(67) 



where we made subtraction at zero renormalization scale C = 0- For purpose of 
brevity we have deflned n{y) = nji{y = p\; C = 0) + m^(0), variables x, y, z are 
squares of Euclidean momenta such that z = x -\- y — 2yjxy cos9 

In our numerical treatment we put m^(0) = 1 in arbitrary units. The coupling 
strength Cg is then scaled through the dimension-full coupling constant g which 
is taken in the units of m(0). 
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Figure 2. The renormalized mass plus selfenergy for the coupling strength Cg — = 
0.5m^(0). Minkowski result is labeled by Unitary Equations, EU stands for the solution of 
Euclidean SDE, PT means one loop perturbation theory. 



Within the same choice of renormaUzation scale C and the same choice of 
renormahzed mass the Minkowski problem has been solved to. The resulting 
'total selfenergy' function i7R(p^;0) +m^(0) is plotted in Fig. I^lfor the coupling 
strength Cg = 0.5m^(0). The Unitary Equations as well as the Eq. (|67|) has 
been solved by the iterations without some peculiar troubles. We use Gaussian 
quadrature method for the integration. Note only that in order to calculate n{zi) 
in (|H7|) we need to interpolate (extrapolate) between the 'fitted' points n{xi). To 
achieve good numerical accuracy in the Minkowski approach we have found it is 
necessary to use relatively large number of mesh points. This is a consequence of 
the presence of Heaviside step functions in the kernel and we take 800 of mesh 
points in the case of figure's calculation. The value of obtained physical mass is 
w? = 0.92m^(0) for Cg = 0.5m^(0). Furthermore, we compare with the one loop 
perturbation theory. These results are calculated with the physical pole mass 
as it has been determined from Unitary Equations. Note that determining pole 
mass from Eq. (|(i4|) in a framework of perturbation theory would lead to the 
small negligible shift of pole mass and threshold as well. 

3.2 The relativistic bound-states 

In quantum field theory the two body bound state is described by the bound 
state vertex function or, equivalently, by Bethe-Salpeter amplitudes, both of 
them are solutions of the corresponding (see Fig. IS} covariant four-dimensional 
Bethe-Salpeter equation [HHI- In the so-called ladder approximation the scattering 
matrix is given by the sum of the generated ladders. When the more realistic 
models are considered (hadronic physics for example) one is forced to use dressed 
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correlation functions instead of their bare approximations. 

Here we are considering simple super-renormalizable model with cubic in- 
teraction (/>?^3 ;i = 1,2 of three massive fields ( thus this model can be called 
massive Wick-Cutkosky model), where in the approximation that allows us to 
exploit the results of previous Section. Having the propagators calculated within 
the bare vertex truncation of the SDEs we combine them with the dressed ladder 
approximation for the scalar s-wave bound state amplitudes. To do this explicitly 
we follow the treatment described in ^] where the spectral technique was used 
to obtain the accurate results directly in Minkowski space. In the mentioned 
paper the analytic formula has been derived for the resulting equation which 
significantly simplify the numerical treatment. For the extensive but certainly 
not exhaustive review of up to date known or used methods in scalar models see 
also 

For our purpose we utilize the PTIR already reviewed in the Section 3. Very 
similar to the SDEs treatment reviewed in the previous Sections, here the Bethe- 
Salpeter equation written in momentum space is converted into a real integral 
equation for a real weight function. This then allows us to treat the ladder Bethe- 
Salpeter equation in which all propagators (of constituents and of the exchanged 
particle) are fully dressed. This is achieved by the implementation of SR of the 
propagator which is determined by the appropriate set of SDEs. Having solved 
equations for spectral functions, one can easily determine the Bethe-Salpeter 
amplitudes in an arbitrary reference frame. 



3.3 Solution of Bethe-Salpeter equation in Minkowski space 

The Bethe-Salpeter amplitude for bound state {4'i,(p2) in momentum space is 
defined through the Fourier transform of 

{0\TM^^)Mx2)\P) = e-^^-^l ^e-^-c|>(p, P) , (68) 

where X = rjixi + 1J2X2 and x = xi — X2, so that xi = X + i]2X, X2 = X — rjix. 
Here pi^2 are the four-momenta of particles corresponding to the fields (j)i^2 that 
constitute the bound state {4>i,(j)2)- The total and relative momenta are then 
given as P = pi + p2 and p = {rj2Pi — ??iP2)i respectively, and = M^, where 
M is the mass of the bound state. Finally, P ■ X + p ■ x = pi ■ xi + p2 ■ X2- From 
now on we will put ryi = r]2 = 1/2 for simplicity. 

Introducing the scalar Bethe-Salpeter vertex function F = iG^^G2^^, the 
homogeneous Bethe-Salpeter equation for a s-wave bound state reads 

Vip, k; P)Giik + P/2)G2{-k + P/2)r{k, P) . (69) 

The bound states appear as poles of the scattering matrix at total momenta 
< < (mi + 7712)'^ from which the constrain on the normalization condition 
follows (for the derivation see for instance [SSI)- 

In the original Wick-Cutkosky model [HZI the exchanged boson is massless in 
the kernel V and no radiative corrections are considered. This model is partic- 
ularly interesting because it is the only example which is solvable exactly PT|. 
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p+ P/2 





-p+P/2 

Figure 3. Diagrammatic representation of the Bethe-Salpeter equation for the bound state 
vertex function. 



For this model the s-wave bound state the PTIR reduces to the one dimensional 
expression: 

1 

rip, P)= I dz- 



-1 



— {p"^ + zp ■ P + p2/4) — ie 



where the function p{z) satisfies the equation 



(70) 



dzV^^\z\z)p{z), 



(71) 



^ e{s{z-z'))Ts 
^ 2{m? - S) 

s=± ^ ' 



S 



l-z 



12 



l±z' 
l±z' 



The resulting weight functions p{z) for various fraction of binding r/ = v P"^ /2m 
are displayed in Fig. |1J 

When the mass of exchanged particle is nonzero or even when the other 
diagram is considered in the kernel of Bethe-Salpeter equation then the appro- 
priate PTIR for the Bethe-Salpeter vertex function is two dimensional [57| (see 
alsoJSSl)) 

1 oo 

pW(a,z) 



r{p,p) = J dz 

-1 



da 



(72) 



[a - [p"^ + zp- P + PV4) - ie]'' ' 

Thus the equation for original Wick-Cutkosky model is regarded as only the 
special case of ()72() where the weight factorizes and is singular in a. However 
here, all the particles are massive thus the positive integer n represents a dummy 
parameter constrained only by convergence of the solution. The choice n = 2 was 
found to be a reasonable one. 

In this case the Bethe-Salpeter equation can be converted to the following 
real integral equation for the real spectral function: 

1 oo 



9 



dz 



daV^^\a,z';a,z)p^^\a,z), (73) 
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The interaction kernel is V{p, k; P) = g^Gs{p — k) in the dressed ladder approx- 
imation, where all the propagators are understood to be fully dressed. For this 
case the kernel in spectral Eq. H73() yt^l (a', z'; a, z) has been derived in The 
explicit formula reads 



3 

T 7=1 i=± 



e{xi)e{i - xi)e{D) 

2J(a,z)2 



TJ{a, z) 



{l-Xi)\E{x^,S,a') 



sgn [E{xi, S, a')) ln(l - Xi] 



(74) 
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where the functions entering the above relation are 

a' - as - i?T ± 

D = {Rt - a' + as)"^ - iasiRr - S) , 
E(x±,S,a') = a ^ — S, 

«! + a2 ai — 02 / 
Rt = J{a,z}T^ \ z , 

01+02 01 -02 

J{a,z) = a z, (75) 

and where the symbol J2t should be understood as 

E /(^) = /(O) - - ^')/(r+) - eiz' - z)f{T^) , (76) 

T 

where / is an arbitrary function. 

The spectral function cTj in Eq. 1)741) corresponds with the SR of i-particle 
propagator of i-field in the considered and relevant part of the Lagrangian 

^{strong) = - 9(t>t Ml, " | '^'203 • (77) 

These were obtained from the Unitary Equations. Note that n = 2 in DRs (|57l) 
because of two subtractions were made, i.e. 



iT,R(/) = 77i(/) _ 77,(0 - (/ - CO ^^^^TT^\p^=Q , (78) 



since the on-shell renormalization prescription Q = mf with unit residuum was 
chosen. 

The derivation of DRs for all three 77j/j is very straightforward and follows 
the lines of previous Section (for some details see [SB)- The absorptive part of 
n is given by pi/vr, these are 

f 

Pi{<^) = J^^:^ J dad/3X{a,/3;uj)a3{a)ai{/3) , i = l,2, 

P3{^) = ^ dadl3X{a,(3;uj)ai{a)ai{(3) , (79) 

where X is defined by Eq. (j^ . 

In the ladder approximation the obtained spectrum agree with the one ob- 
tained by other techniques [HHl EOl IMl 1^ • In dressed ladder approximation there 
is no other reliable result published in literature to be compared with. However 
the qualitative agreement with the paper jSH] was found: The critical value of 
coupling constant gc (defined by the constraint for the renormalization constant: 
Z{gc) = 0) as it has read from SDEs gives the physical domain of applicability 
of Bethe-Salpeter equation. 
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It is interesting to note here that within the use of Feynman-Schwinger rep- 
resentation method there was found in [IJ that the crossed boxes are dominates 
in scalar models studied here. For a comparison of different technique, including 
Bethe-Salpeter see the result in this paper. The observed dumping of bound state 
energy is due to the bose statistic of the matter field in the model considered here 
(while there is an expected cancellation between higher order skeletons when one 
consider the bound states of fermions). 

At this place we should mention the so called 'ghost' solution of Bethe- 
Salpeter equation. These, probably unphysical solutions (see [SHI and references 
therein), have been found only for rather strong coupling in ladder approxi- 
mation. However, the couplings below the critical one allow only solutions for 
relatively weakly bound states and neither 'ghost' solution of the Bethe-Salpeter 
equation was found. It suggests the domain of validity of SDEs restricts the solu- 
tions of the Bethe-Salpeter equation in such way that only the normal solutions 
exist. 

3.4 Elastic Electromagnetic Form Factor 

Although the elastic form factor represents a simple dynamical observable, its 
Minkowski calculation represents a nontrivial task. For this purpose we consider 
the massive Wick-Cutkosky model given by a Lagrangian gauged as follows: 



where the covariant derivative is = — ieAfj_. 

The electromagnetic form factors parametrize the response of bound systems 
to external electromagnetic field. The calculation of these observables within the 
Bethe-Salpeter framework proceeds along the Mandelstam's formalism |42j. In 
the form factor calculations the effects of scalar dressing were not taken into 
account. 

The current conservation implies the parametrization of the current matrix 
element in terms of the single real form factor G{Q'^) 



where = —q'^, so that is positive for elastic kinematics. 

The matrix element of the current in relativistic impulse approximation is 
depicted in Fig. |21 The matrix element is given in terms of the Bethe-Salpeter 
vertex functions as 




(80) 




(81) 




(82) 



where we denote P = Pi and represents one-body current for particle (pi, which 
for the bare particle reads jiiPf,Pi) = +Pi^ where Pi,Pf is initial and final 
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f\ p+q+P/2 p+P/2 A 



P+q 




-p+P/2 



Figure 5. Diagrammatic representation of the electromagnetic current bound state matrix 
element. 



momentum of charged particle inside the loop in Fig. El i.e., pi = p + P/2,pf = 
q+p + P/2. 



The main result, as far as charge form factor is concerned, is the rewriting of 
the rhs of Eq. H82|) directly in terms of the spectral weights of the bound state 
vertex function. It allows the evaluation of the form factor by calculating the 
dispersion relation: 



without having to reconstruct the vertex functions r{p, P) from their spectral 
representation. Unfortunately the derived expression for pc in the DR H83|) con- 
tains some additional integrations which the authors were not able to remove. 
Due to this and in order to avoid sizable numerical errors the results have been 
calculated only for positive . However a dramatic and exciting changes can be 
expected at timelike momentum regime. This remained unexplored until now. 

4 Strong coupling QED 

In this Section we investigate SDEs in 3+1 dimensional QED. We follow the 
prospects presented in the papers |13[ I14j , in which the solution for this theory 
in Euclidean and Minkowski space were compared for the first time. 

The first part of this Section is devoted to the solution for the electron propa- 
gator in ladder approximation, in the second part an extension to the unquenched 
case (with the photon polarization included) is made. In the latter case the run- 
ning coupling is considered self-consistently: the fermion mass and photon polar- 
ization function have been calculated by solving the coupled SDEs for electron 
and photon propagators in Landau gauge. In this paper we do not use the vertex 
ansatze as in the paper |14j but we consider the first term of skeleton expansion 
(i.e. bare vertex is used). Having small quantitative effect, because of Landau 
gauge, this not only drastically simplified the calculation effort but allows us to 
perform many integrations analytically. Likewise in the Section where we discuss 
a scalar model we write down the Unitary equations where the principal value 
integrations are excluded. This is a clear advantage when compared with earlier 
version of Unitary Equations (T^ITl] where this possibility was overlooked. 




(83) 
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In both approximations we are looking for solutions with zero and nonzero 
electron mass in Lagrangian. In the first case the chirally symmetric solution 
(with massless electrons) always exists. In the Euclidean formalism we obtained 
agreement with many previous studies |431l44lll5] . for sufficiently large coupling a 
also nontrivial solution for the mass function. On the other hand, no such solution 
was found in our spectral Minkowski approach. This is a strong indication that 
the dynamical fermion mass function has in this case a complicated (complex) 
structure, which is not reflected by assumptions of our spectral Ansatze. 

We will be dealing mostly with the strong coupling regime, which is far from 
"real- life QED", for which uqed (in experimentally accessible energy region) is 
small and use of perturbation theory is fully justified. Of course, all our solutions 
agree in weak coupling limit with the perturbation theory. 

As it is already stated in the introduction the basic motivation is to decode 
the timelike structure of nonperturbatively obtained Green's functions. For this 
purpose the models considered here are an ideal playground. However, in this 
Section the reader can find at least two other reasons to consider strong-coupling 
QED dynamics. 

The strong-coupling Abelian dynamics was considered as one of the candi- 
dates for explanation of electroweak symmetry breaking |46,. Such or a similar 
models can approximate the dynamics in the strong coupling sector in various 
Technicolors model, e.g.. Slowly Walking Technicolors [471 1481 l5Uj . for which 
ladder approximation of QED-like theory seems to provide a test background for 
modeling of the (Techni)lepton propagators (for a possible indirect experimen- 
tal evidence of Technicolor-like model see [S^)- Moreover, the strong-coupling 
Abelian model is not only a suitable playground for studies of supercritical phe- 
nomena like dynamical mass generation, but it is suitable also for investigations 
of the analytical structure of the fermion propagator. It was 30 years ago when 
Fukuda and Kugo |H2] observed the disappearance of the real pole of the fermion 
propagator in the ladder QED and it was argued that this is a signal for confine- 
ment of the fermion. Again, the main motivation is the simplicity of the model, in 
alternative more complicated models the nonperturbative phenomena are even 
more difficult to study quantitatively. We believe that it is worth reviewing the 
recent progress in strong coupling QED. 

Our attention is mainly focused on the SDE for the fermion propagator, 
because no nontrivial nonperturbative effects (i.e. effects not known from per- 
turbation theory) were found in the photon propagator. However, let us remind 
the solution of perturbation theory renormalization group equations for the run- 
ning coupling [53! ) which in two loop approximations exhibits Landau singularity 
at the scale 




[21 



(84) 



where /3j are the coefficients of QED beta function 




+ ... = /3iefj + /32e| + ... 



(85) 
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giving us numerically vIl — 10 GeV if only the electrons are considered. For the 
particle spectrum of the standard model one gets vIl — lO^^GeV which is many 
orders above the Planck scale and hence probably without any physical relevance. 
In the minimal Supersymmetric Standard Model with four Higgs, which offers 
a solution to the strong CP problem, the Landau pole moves down to Ai ~ 
10^'^ GeV. Thus the Landau pole is by no means academic and to resolve this 
problem one has to consider the nonperturbative formulation of QED. For this 
purpose we will consider coupled SDEs in unquenched QED but here rather in 
the strong coupling regime where the Landau singularity is expected at a much 
lower scale. 



4.1 QED SDEs 

The SDEs for one particle irreducible Green functions can be derived by applying 
the functional differentiation on effective action /^[-^c, "^c, ^c], where '0c(^c) is 
semiclassical vacuum expectation values of fermion (photon) field in the presence 
of the fictitious sources. For instance, the inverse of photon propagator is obtained 
as: 

where all effective fields are set to zero at the end of the derivation. The details of 
derivation of SDEs can be found in the in standard textbooks, (see e.g. pp. 289- 
294 of Ref. [SS]) but also |^|S5)- Here we review the SDEs for propagators in 
momentum space for completeness here. The SDEs for the inverse of photon 
propagator G^^'^ and for the inverse of electron propagator S read 



(G-i),.(x, y) = -^L±i;i£l|^ , (86) 



+ nf^'iq), (87) 



n^''{q)='^J^^y^=ie^Tr J S{k)r {k,k - q)S{k - q) , 



S-\p)=Pf,j^-mo-U{p), (89) 
m=-^^^^=ie' I ■^rS{k)r'^{k,p)G.,{k-p), (90) 

where n'^'^{q) is the vacuum polarization tensor, IJ{p) is the fermion self-energy 
and the blobs in the diagrams represent photon-fermion vertex function satis- 
fying its own SDE which can be sketched as in the Fig.lHl Further, the parameter 

in Eq. ()87() follows from the common choice of covariant gauge fixing term 
-(1/2^^) {df.At'f of the QED action g4j. 

From the conservation of the gauge current the following Ward-Takahashi 
identity for the photon-fermion vertex follows: 

ip-i)^r^{p,i) = s-Hp)-s-\i), (91) 
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Figure 6. Diagrammatical representation of the Schwinger-Dyson equations for vertex function. 
Like in the diagrams for selfenergy the blobs (box) represent the exact vertex and the internal 
lines represent the exact propagators- wavy line stands for the photon and the straight line 
labels fermion propagator. 



while the same imphes transversahty of the polarization tensor: 

g,i7^''((Z) = 0, (92) 

noting here that Eq. (|92|) also follows from Eq. (|88j) using Eq. (|91j) provided the 
divergent integrals in (|88() are regularized by a translation invariant manner. 

Therefore, in this case, the vacuum polarization tensor can be parameterized 
in terms of single scalar function 7T(q^): 



n'"'{q)=q'^ 



After the inversion, it yields for the photon propagator 



where we define the photon renormalization function H{q^) as: 
From H94() one finds that the bare photon propagator is given by: 



The general form for fermion propagator ()89() is 



or alternatively: 



S{p) 



F(p2 



M(p2 



(93) 



(94) 



(95) 



(96) 



(97) 



(98) 



where F{p^) = A~'^{p^) is called the fermion wavefunction renormalization and 
M{p'^) = B/A is the fermion mass function. Clearly the fermion propagator for 
a free fermion field or bare fermion propagator is given by i7 = in Eq. ()89|) . 
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4.2 Physical constraints on the solution 

The model discussed in Section |31 was too simple and thus does not render some 
very important issues encountered when one deals with more realistic models, 
e.g. the models for which we use a gauge theory concept as a starting point. 
Therefore we discus here the physical criteria which make the solution of SDEs 
(and any other nonperturbatively obtained solution of Green's functions) physi- 
cally meaningful. 

In Quantum Field Theory an experimental observable is given by the tran- 
sition probability associated with the appropriate part of the S-matrix. When 
dealing with a renormalizable gauge theory the S-matrix element describing a 
given process should be completely independent on the renormalization scale (it 
is renormalization group invariant) and on the particular choice of gauge-fixing 
parameters (it is gauge-fixing independent). Furthermore - as shown above - the 
GFs should satisfy the Ward identities. For a sufficiently small coupling constant 
the usual perturbation theory in conventional gauges offers a safe way to satisfy 
all the required invar iance. It is well known that the gauge variant (unphysi- 
cal) parts usually involved in the original GFs cancel out in a given S-matrix 
element. To this end, there is recent progress in the so-called Pinch Technique 
j55| EE] which offers gauge-fixing independent Green's functions to all orders of 
perturbation theory both for the Abelian and nonAbelian case. It is also re- 
markable that similar to the Background- Field Method the Pinch Technique GFs 
automatically satisfy the naive Ward identities rather than more complicated 
Slavnov- Taylor identities. 

The problem of gauge invariance in the nonperturbative SDEs formalism 
is not yet fully solved. Adopted techniques should maximally respect the gauge 
identities such as Ward identities, Nielsen identities [57j and Landau-Khalatnikov 
transformation |58| I59j when class of linear covariant gauges in QED is used. 
Clearly, the best we could have from this point of view is the systematic method 
of building the renormalized blocks of GFs which are already free of unphysical 
information. Very promising but rather complicated method is to circumvent the 
gauge fixing problem by the use of stochastic quantization. This unconventional 
idea is based on the introduction of a fifth stochastic-time coordinate, while the 
usual observables are obtained in the equilibrium. (The stochastic quantization 
and its relation to the usual methods is up to date exhaustively reviewed in 
j6Uj). Since contrary to the effective action approach here, this method does not 
require gauge fixing at all and hence it offers also the gauge fixing invariant 
GFs. Further, it is important to note that within the certain approximations 
the stochastic quantization yields the same Yang-Mills SDEs as conventional 
approach in Landau gauge [HSl- The second promising method could be above 
mentioned Pinch Technique which was originally developed in the context of 
SDEs |Snin21 in 1980's. However, in this strategy the requirement of multiplica- 
tive renormalizability has been ignored and to the best of author's knowledge 
from that time no progress has been made in utilization of Pinch technique in 
the SDEs formalism. 

The gauge technique [SU ES]) being based on the spectral representation 
offers GFs satisfying Ward identity in the whole Minkowski space. Within the 
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context of gauge technique R. Delbourgo in Ref. [OS] highlights some of the 
related problems hindering the studies of SDEs. Among these, the first one is how 
well the Green's function solutions obey the Landau-Khalatnikov transformation 
and the second is how closely do the nonperturbative answers coincide with the 
perturbation series if we expand them in the power of couplings. 

Making a simple Ansatz for a gauge vertices one can ensure Ward identity 
irrespective of analytical properties of GF's. In QED the simplest solution is 
the Ball-Chiu vertex [67^, while any improved Ansatz can differ at most in the 
transversal part of vertex function. For an exhaustive list of the modeling vertices 
in QED see the papers |681 169j . However stressed here that the recent numerical 
study [ZOI of the gauge dependence of the chiral condensate leads us to the 
conclusion that the Ward identity alone is not sufficient to ensure the gauge 
independence of the physical observable. Further, it is argued that it is essential 
to apply full Landau-Khalatnikov transformation to the dynamically generated 
mass function as also advocated in [TJ. This step seems to be necessary in order 
to obtain a reliable results in other gauges then Landau gauge. 

With the exception of the gauge technique to solve the system of SDEs within 
the help of spectral technique means to calculate complicated formula which 
should lead to the evaluation of dispersion relation for proper GFs at the end. 
Therefore we instead considering some Ansatz for the vertex function we rather 
use the (two simplest) skeleton expansion of self-energies and our calculation is 
performed in conventional Landau gauge where = 0. Such skeleton expansion 
can be systematically improved by considering further and further terms in the 
multiloop expansion. In agreement with the issues made in the paper j72] we 
assume that the usage of Landau gauge should lead to the approximately gauge 
invariant GFs and gauge-fixing independent observables. 

The renormalizability means that divergences may be removed (order-by or- 
der in perturbation theory) upon the field redefinitions 

^ ^ yZsV', ^ v%Af„ (100) 

and vertex renormalization 

7^ ^ Zi7^. (101) 
If the renormalization scheme respects the Ward-Takahashi identity, i.e. 

Zi = Z2 (102) 

in the exact theory we say that QED is multiplicatively renormalizable. 

To summarize, we impose the physical criteria, which will allow one to deter- 
mine the solution of SDEs to be physically meaningful (in general renormalizable 
gauge theory). 

• The theory should be multiplicatively renormalizable, which implies: 
-the number of possible subtractions for a given GF (in momentum sub- 
traction scheme). 

-the relation between various renormalization constants (and hence the co- 
efficient of subtracting polynomials). 
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-renormalized GFs satisfy WTIs. 

-various possible renormalization schemes are physically equivalent 
-the position of the pole in the propagators is scheme-invariant. 

In addition, for the theory to be tractable in our Minkowski approach the 
following is assumed or required 

• The propagators are analytical functions in the whole complex plane of the 
momenta except the real positive axis, which implies: 

- the vertices satisfy integral representation known from perturbation the- 
ory (here it is assured by skeleton expansion of GFs) 

- the proper Green's functions satisfy dispersion relations 

However, note here: neither the existence of a pole nor the positivity of the SR 
is needed. 

4.3 Ladder approximation of electron SDE with nonzero bare mass 

In ladder approximation the vacuum polarization is neglected: the photon propa- 
gator is taken in its bare form and the photon-fermion-fermion vertex is approx- 
imated by the pure 7^ matrix. The relative technical simplicity of calculations in 
the ladder approximation for electron SDE is an attractive feature rendering very 
interesting study without some large computational effort. It allows transparent 
comparison with the Euclidean and the perturbation theory results as well. 

Using bare vertex and bare photon propagator leads to the following ladder 
fermion SDE: 

S-\p) = Z2[^ - mo] - U{p) , 

m = f ^,G^freeik)l,S{p-k)^.. (103) 

Further we introduce the mass renormalization constant and/or the mass 
counterterm 5m 

mo = Zm m{C) ; (5m = mo - m(C) , (104) 

relating the renormalized mass m{C,) at the renormalization scale C, to the bare 
mass mo and recall that Zi = Z2 = 1 in ladder approximation in Landau gauge. 

4.4 Euclidean formulation 

First of all, in order to make a careful and constructive comparison with the 
standard Euclidean results we review the basics of the standard Euclidean for- 
mulation j52| 143 1 144j results presented in the literature. Taking the traces Tr 
and Tr in the first Eq. H103() the SDE is transformed into one equation for the 
mass Me = Be ■ 

Ae{x) = 1, 

Be{x) = mo + ^ / dyK{x,y) f^jj^) (105) 


2y 



V. Sauli 



31 




Figure 7. Spacelike solutions for mass function of the electron propagator as have been obtained 
within spectral (S) and Euclidean (E) formalism. 



where Wick rotation and angle integration have been done and where we have 
used Z2 = 1, X = p'^ = —p'^, y = k'^ = —k'^ . The renormahzed Eq. (|1U5|) then 
reads 

00 

Be{C,x) = m(C) + ^ / dyV{Cx,y) f^J'J^ (106) 
4vr7 y + B^iCy) 

V{C,x,y) = K{x,y)-K{C,y). 

Recall here that Eq. (|l()5j) with zero bare mass (mo = 0) provides nontrivial 
solution only for a > etc (tlie so called supercritical QED, Qc is critical coupling) 
while BijP') = for a < Uc (chiral symmetric phase). Recall also here that the 
supercritical solution requires the implementation of a UV regulator Ah owinig 
to Miransky scaling [JSl 1^ : 



B{Q)^AhJi--^ (107) 
V a 

where Ah is the naive (translation invariance violating) integral momentum cut- 
off. 

In this paper we solved the SDE H106|) with the renormalization choice M{C, = 
-10^) = 10. The results are displayed in Fig. [3 It is not surprising that for 
a > ac the obtained results do not agree with perturbation theory at all. In 
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agreement with the studies j52[ I44j the expected damping of the mass function 
to its negative values is observed in supercritical phase of QED. 



4.5 Minkowski formulation 

As in the scalar theory discussed earlier, we use the SR for the propagator and 
arrive at DRs for selfenergy. Although the results of this Section are presented for 
ladder approximation QED the renormalization wave function A will be formally 
kept during the calculation for a possible more general future studies. There are 
two independent spectral functions in SR for QED fermion propagator: 

S{p)= [ds^^^l±^ = -^+ [ds^^4^1±^, (108) 
J p-^ — s + le p — m J p"^ — s + le 

where m is the physical electron mass. As usual the propagator is split into a real 
pole part and a part that is generated by the interactions. To avoid confusion 
from the beginning, it is stressed here that the assumption ()lU8p does not imply 
(but does not exclude) real or complex pure pole singularity structure of the 
electron propagator. Thanks to the masslessness of the photon, the branch point 
corresponds to the physical mass of the electron. The gauge-dependent complex 
singularity appears after the limit p'^ — > m^, while under the threshold p'^ 
the propagator is assumed to be a real function in any case. This known fact is 
reflected by the threshold behavior of the continuous spectral function a^^s (in 
the previous scalar case the spectral function was zero at the threshold, however 
here ay^s{fn-) / in general). 

We split the unrenormalized self-energy H103|) into its Dirac vector and Dirac 
scalar part 

E{p) =^a{p^) + b{p^) . (109) 

The self-energy in Eq. ()1U3() is only logarithmically divergent and one subtrac- 
tion: 

aRiC-.p") = a{p') - a(C) = / dco . /-^^y -^ , (no) 

bniCp') = Hp') - 6(C) = / du /-^"^.T^^ . (Ill) 

J (p^ - uj + ie){uj - Q 

is sufficient to renormalize the self-energy: I!r{C,,p) =^aR{(;p'^) + 6_r(C;P^)! 
where two scalar functions a, b satisfy the DRs written also above. 

Relatively straightforward calculations OSl lead to the results for the 
imaginary parts of the functions (|110() . (|111() with the result p„ = while the 
weight function ps reads 



= -3 (£) 



rm (1 - ^ ) + / daas(a) ( 1 - ) 
V P J Jm2 V P J 



©(/ - m^) . 



(112) 

As we already saw in the previous Section the peculiarity of the ladder ap- 
proximation electron SDE in Landau gauge gives for the renormalization wave 
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function A{(^;p'^) = 1 + aji{(^;p'^) = 1 and renormalization group invariant mass 
function is 

M(p2) = B{p^) = m(C) + bniC; p") . (113) 
The Unitary equations can be most easily derived from the fohowing identity 

S-^S = l, (114) 

where we use the SR (jlOHf) to express fermion propagator S and SDE in DR 
form for the inverse, which reads 

=^ - m(C) - / ds ^^ , ^^v;j.-c 
J p"^ — s + le 

For this purpose we introduce the fohowing convenient notation: 

g = ReS = — \- P. ds^ ; (t(s) =^fTt,(s) + cJs(s) , 

p — m J p'^ — s 

d = ReS'^ =^ - m(C) -P. [ ds-^^^ ; p{s) =^p,{s) + ps{s) , 

J p s s (, 

(116) 

which ahows us to rewrite the identity ()114() as 

ig-iTTa{p^)){d + i7rp{p^)) = l, (117) 

where we have used the well known functional identity for distributions (|59|) . 
Comparing the real and imaginary parts of (|117|1 we arrive at two coupled equa- 
tions: 

gd = I- 7r^(T(p^)7rp(p^) , 
<T(/)d = gp{p^). (118) 

Still in the matrix formulation introduced in Eq. ()116j) we can write 

i.e. g is eliminated. Owing to this fact the principal value integrations over the 
unknown & are not presented from now on. The Eq. H119|) represents matrix form 
of Unitary Equations. For the purpose of completeness we will write them down 
also in terms of cr^^s explicitly. To this end we define two scalar functions Ct,(p^) 
and Csip^) such that 

d2 + 7rV(P^) =:^c„ + c,, (120) 

which implies 

Cs=p^A^ +B + t:'^{p^pI + pI) ; = -2 AB + 2tt'^ p^p^ , (121) 
where we used shorthand notation for the real parts of the functions A, B 



A = Re^(/) = 1-P. j ds 



p^ — s s — C ' 
p^ — s s — C 



B = ReB{p'^)=m{C) + P- I ds-^^^-^. (122) 
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Then "projecting" Eq. by Tr and Tr ^ we get, after a little algebra, the 

desired Unitary Equations: 



a,, 



CsPs 



Cy 



^vPs ^sPv 
p2 d - cl 



(123) 



which are nonzero above the threshold, i.e. for > rm?. As in the previously 



discussed scalar model, the form of the Unitary Equations (|123j) does not de- 
pend on the details of interaction but on an analytical assumption - a spectral 
representation for the fermion propagator. Clearly, the Unitary Equations are 
formally equivalent to the Unitary Equations presented in the Ref . PS| j but here 
with the main difference that all principal value integrals to be solved numeri- 
cally are eliminated from now on. In our approximations the function A = 1 and 
the real part of the mass function above the threshold is 



m(C) + 



4tt 



r m 



1 



1 



+ 



daaJa) 



4lln 



m 



a 



In 



p^ 



1 



(124) 



As we mentioned, the threshold singularity need not be real valued and the 
convential procedure of extracting a real number r by on-shell differentiation of 
the self-energy function is inconsistent. However, the solution of this puzzle is 
very simple, recall here that the fixed renormalization uniquely determine the 
whole propagator function and thus the real pole part residuum r in p()8|) as 
well. Putting for instance p = in (|114l) we get particularly simple formula for 
the residuum r 



m 



m{0 + L 



m? x{x—C,) 



dx 



a^[x m 



-dx . 



X 



(125) 



The physical electron mass is defined by Eq. ReS ^(m) = or equivalently 
M{m) = m. Using the DRs (|110|) . (|111|) the desired relation reads: 



m = m{() + 



as- 







{m? — s){s — C) 



(126) 



However, here Pv^sifn?) = is necessary in order to have the correct definition 
of the physical mass. (This statement is clearly fulfilled in perturbation theory, 
see the first term in Eqs. (|112|1 . (112411 for this purpose.) 

The Unitary Equations ()112I) . H123() have been solved by the method of it- 
erations. The residuum (|125|) . physical mass (|126j) and the mass function (|124|) 
have been evaluated at each step of the iteration procedure and substituted in 
the Unitary Equations until convergence is achieved. 

To compare with the Euclidean space SDEs solutions we use the same renor- 
malization condition in the both approaches, i.e. -B(— 100) = 10. Having all 
the numerical solutions stable and making comparison between Euclidean and 
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Minkowski results we have found that they perfectly agree when the coupling is 
small enough. The Unitary Equations resulting spacelike solutions for the mass 
function are compared in the Fig. [7J In Fig. |H1 the Unitary Equations solutions 
are given for timelike momenta. The meaning of the origin of small numerical 
difference for strong couplings is not satisfactory understood. We can speculate 
here that the analytical spectral representation Ansatz does not fully capture 
the structure of the fermion propagator but we hope this point will be clarified 
in the future. 

4.6 Approximate solution of dynamical chiral symmetry breaking in Minkowski 
space 

In explicit chiral symmetry breaking case, in which the nonzero electron mass 
exists from the very beginning, both approaches- Minkowski and Euclidean- 
offer approximately the same results in the coupling constant domain where 
both solutions were obtained (in the spacelike momentum region, of course). 
These solutions were reviewed and compared in the previous Sections. 

Until now, we have avoided the phenomenon of spontaneous chiral symme- 
try breaking which plays an important role in particle physics. The numerical 
methods used for rainbow Landau gauge approximation in the paper enable 
to find a non-trivial solution of Unitary Equations for zero fermion Lagrangian 
mass. The numerical solution tended to fail near critical value ac = vr/3 in rain- 
bow approximation. In that case, approaching the coupling to its critical value 
the Unitary equations become less numerically stable. It was already observed 
that better stability can be achieved when nonzero mass is added into the gauge- 
boson propagator. The authors of the paper [75] were looking for the Minkowski 



36 



Implications of analyticity to. 



solution of the fermion SDE i7 = J KS with the effective kernel chosen as 



Kiq,p)ijki 



a P 



a/3 , 
i/ ~r 9 



1 



— rrij^ + ie 



(127) 



In the case when uib « tha the mass m b can be understood as an infrared reg- 
ulator while my[ represents ultraviolet (Pauli-Villars) one. For this simple model 
the gap equation has been solved in the Euclidean space and in Minkowski space 
as well. The nontrivial chiral symmetry breaking solutions have been compared 
for the first time. Within the spectral function Ansatz H108|) the strategy of 
Unitary equations (|123|) has been used for the later purpose. In ^ = 1 approx- 
imation the explicit evaluation of finite selfenergy leads to the following result 
for the absorptive and dispersive part of the mass function 



a 



(47r) 



dxasix) [Xo{u;;m%,x) - Xo{u;;m\,x)] , 



Mip^) = ^ j dxasix) [J{p\x,ml)- J{p\x,m\)] , (128) 

where we have labeled 

0(-Ap)v^ 



0(Ap)7(A, 

p2 



TT p — y — z 
- + arctg ^— 



+ ^ln(16yz) 



In 



Ar, 



p^ -y-z + 5^ 



P^ - y - z + JXp 



0(AoVA^. 
H K m 



-y- z + Vx^ 



(129) 



and the following abbreviations is introduced: 

Ap = X{p'^,y,z), 

\/Ao = \y-z\, 
T = iVy + V^? 



(130) 



The function X is defined in ((24 

These equations together with ()123() comprise the coupled set of integral 
equation that have been solve numerically in jJH]- Not strangely, the critical cou- 
pling Oc ~ vr/3 when ms « tha which is expected from the study of rainbow 
QED with hard cutoff. More interestingly, the analytical continuation of (spec- 
tral) timelike solution to the spacelike regime is in rather good agreement with 
the Euclidean one, especially when the coupling is close enough to the critical 
value Oc- Furthermore, the trivial solution was found bellow Oc in both treat- 
ments. However, there is apparent deviation when the coupling increases. This, 
together with the authour's observation of the deviation from assumed analyt- 
icity implies that spectral solution is not the exact and true solution of original 
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SDE. It is suggested that deviation of 'Minkowski spectral' solutions from cor- 
responding Euclidean ones reveals that the employed Ansatz, which fixes the 
analytical structure of propagator in Minkowski space, is not adequate for the 
case of the dynamical mass generation. It seems likely that with coupling increas- 
ing above its critical value, the pole of the propagator, located for subcritical case 
on the real axis, moves gradually into complex plane (i.e. receives non-negligible 
imaginary part) and the standard Lehmann representation, employed for our 
calculations, is not valid anymore. 

However, checking the analyticity in the whole Minkowski space and also com- 
paring spectral solution with the Euclidean ones at spacelike, one can conclude 
that the solution of Unitary Equations near the critical coupling is a rather good 
approximation. This approximative character become worse and worse when the 
coupling flows from its critical value, while the agreement is perfect at critical 
coupling (however were the result is trivial). At the end we can mention that 
within the Unitary Equations a given model neither provides the similar picture 
to the behaviour of light (say u,d,s) quark mass in QCD (here, the word 'similiar' 
means the shape and slope of the mass function, or much simply the ratio of the 
value of mass function at two rather distinct scales). The analytical property of 
the solution in such a strong case remains unexplored at least by the presented 
method. For an different analytical guesses of the form of the quark propagator 
in modeled QCD see for instance j76j . 

4.7 Unquenched QED - calculation of the running coupling 

In the previous Section we performed a numerical study of the electron mass 
function in the simple approximation to the SDEs: the bare vertex and pho- 
ton propagators were employed. In what follows we extend these calculations: 
we solve the corresponding SDEs in Landau gauge also for the photon polar- 
ization function. We again compare the solutions for propagators in Minkowski 
and Euclidean space and we have found reasonable agreement between solutions 
obtained in these two technically rather different frameworks. 

The 4D QED is trivial theory, in other words the presence of some regulator 
is required for any nonzero value of the renormalized coupling constant. It is 
convenient to introduce nonperturbative invariant continuum regulator function 
which also render finite the ultraviolet divergences. Such program, completed 
years ago, provides invariant regularization across all quantum field theories (for 
the review see the paper |77)). However in this paper we simply use spectral cutoff 
regulator As for the purpose of Minkowski calculation and the hard momentum 
cutoff Ah for the purpose of Euclidean calculation. Ajj cuts the Euclidean mo- 
mentum integrals while As represents the upper spectrum boundaries in GFs SR 
and DR. They have been put to be equal. As = Ah assuming no significant affect 
on the results. By construction the Minkowski regulator does not violate trans- 
lation invariance of the theory while Euclidean hard cutoff Ah do this. The same 
argument is also valid for the gauge invariance of the theory (assuming gauge in- 
variant truncation of SDEs). The observed discrepancy between the Minkowski 
and the Euclidean solutions, although not so large, is the consequence of the 
regularization scheme dependence in the trivial theory. 
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There are many attempts (ZHllBniEIlESIBnilHllHS] to estimate the form 
and its effect of the transverse part Ft of the full vertex F = Fbc + on the 
behaviour of fermion propagator. The eight independent component of Ft is the 
only unknown part of the fermion and photon SDEs system when solely writ- 
ten. The longitudinal part, named Ball-Chiu vertex Fbc is the minimal Ansatz 
consistent with Ward-Takahashi identity H91() (and of course with perturbation 
theory jZH])- To derive dispersion relation for fermion selfenergy and photon vac- 
uum polarization is complicated task even when Fbc is only considered ( see 
|35j ) . therefore, instead of making some complicated Ansatz for the gauge ver- 
tex we consider the skeleton expansions of the Green functions. The skeleton 
expansion can be obtained from the sophisticated iteration of truncated SDEs 
system corresponding thus with approximated effective action. This so called <P 
derivable approximation (or N-PI) of effective action is used in many branches of 
quantum field theory |HnilH3IHHllHnilS3inilE21inniini- For an intimate relation 
between N-PI effective action and the full SDEs see the discussion in [^Sl- It is 
a systematically improvable method where the improvement is achieved by con- 
sidering further and further terms in the loop expansion of the effective action. 
In the skeleton expansion there are only the classical vertices included in each 
skeleton diagram which make the calculation of DR more tractable. For a first 
attempt made in this paper we calculated the one loop skeleton contribution, the 
resulting gap equation corresponds with the SDEs in the bare vertex approxima- 
tion. The study of SDEs where self-energies are considered in two loop skeleton 
expansion is now under consideration. 

Further in order to reduce the complicated hierarchy of equations we explic- 
itly put A = 1. Within several percentage deviation such approximation can be 
justified in Landau gauge herein [^01 113 EH • 

4.8 Unquenched QED in Euclidean space 

The SDEs are solved in Euclidean space after the Wick rotation /cq ikiE is 
performed. Then the loop integrals should be free of kinematic singularities and 
the GFs are found for positive Euclidean momentum kl^ = k'^ + k"^ + k"^ -\- k\. 

Since we use the hard cutoff Ah in our Euclidean treatment the unpleasant 
quadratic divergence in the polarization function n is generated in this case. Due 
to this we follow usual strategy [3 inSl IHZI and avoid the quadratic divergence by 
the use of Bloch-Pennington projector 



Employing this projection in the SDE for the photon propagator we arrive at 
the following coupled equations: 



nE{q^) 





(131) 



nE{x) 



dy 



y 



I 



dd sin^ e 



2y — 8y cos 



i"^ 9 + 6y^yx cos 9 



y + M|(y) 



2 + M|(z) 
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(132) 

where the variables x, y represent squares of Euchdean momenta, z = x + y — 
2y^yx cos 6. 

After the subtraction 

n§iC;p^) = 77s(p2) - 77s(C') , 

M§{C';p^) = Me{p')-Me{0- (134) 

at some arbitrary C,' and after renormahzation, the Eqs. H132p . (|133p have been 
solved numerically. If needed the change of the renormahzation scale choice is 
performed simply by the utilization of the identity 

i7|(C;p2) = i7|(C';p2)_7l|(C';C). (135) 



4.9 Unquenched QED in Minkowski space 

Likewise in the ladder approximation of QED the SDEs written in momentum 
space are converted to the coupled set of Unitary Equations. The Unitary Equa- 
tions H123() . ()125|) . H126|) for fermion spectral functions is already derived in the 
previous Section. They remain unchanged since they follow from general ana- 
lytic structure of the fermion propagator function. The Unitary Equation for the 
photon spectral function is derived bellow. Before doing this we review the DR 
for the fermion selfenergy in the presence of dressed photon propagator. 

In A = 1 approximation we can identify the electron mass function (|98|) as 

M{p^) = mo+ '^U{p) = mo + 



e-jdajdbj ^7 (p_/)2_, 7^ p _ ^7^136) 

which after the renormahzation (add zero of the form m{Q — m(Q into the Eq. 
(|136() and make the subtraction) leads to the DR 

M{p') = Hd ^-^^ +m(C), (137) 

p.iio) = -3 ( -^Y f dadbas{a)aJb) 0{uj - + Vhf{)-^i^) 

\47r/ J io 

where cr^ is the photon spectral function defined by the SR for photon propagator, 
which in the Landau gauge reads 

G'^'iq) = / db \ — , 

^' J q^-b + ie 



a^{b) = r^5{b) + a^{b) . (139) 
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Since gauge invariance is correctly maintained the photon polarization function 
should posses at most logarithmically superficial divergence. Therefore only one 
subtraction is needed and the DR for renormalized polarization function then 
looks like 



duj 



p^{uj){q^ - C) 
{q^ -uj + ie){uj - C) 



(140) 



where we use (not necessarily) the same renormalization scale as in the case of 
electron propagator. 

In order to arrive to the photon Unitary Equation, the real and imaginary 
part of the identity G~IG^"' = da-y is evaluated. This, after using the definition 
of photon propagator (|94() and the SR assumption p39j) leads to the following 
identity: 



+ 



ds a.y{s) 



s + ie 



1 



[l-nn{C;q')] 

For the next purpose we use shorthand notation for the function 

oo 



b.,{q') = l-Renti{C;q' 



I -P. ds 



iq^-s){s-C) 



(141) 



(142) 



noting here that P. integration can be performed analytically at least at one 
(skeleton) loop level. Evaluating the real and the imaginary part of Eq. H141|l 
gives 



P. / ds 



cr^(s) 



q^ — s 



1 -vr 



2 „2 



p^{q^)(T^{q^ 



'■ b^q^) a^q^) - p^{q^) 



ry + q^ P. 



— s 



0, 



(143) 



(144) 



Putting together it yields resulting Unitary Equation for the photon spectral 
function: 



where the functions fj-y and are zero under the threshold q^ = 4m^. 
The residuum of the photon propagator can be calculated from 



(145) 



1-77^(C;0) 



(146) 



In the next text we will sketched the derivation of DR (|14())) . First we briefly 
review the method in its perturbative context. 

In 4 + e dimensions and for spacelike momentum < the one loop polar- 
ization function can be written as 



n{q^) 



4e^ 



3(47r) 



+ — ln{4:TT) + In 



m 

P't'H 



Am? 5 
^ ~ 3 



V. Sauli 



41 



(1 + 2mVg^^ 



4m^ 



-In 



1 + a/1 



4m£ 



1 



4m2 



(^Z3/i,7|, ,(147) 



where ^t'H is t'Hooft scale. The mass-sheh subtraction scheme defines Z3 so that 
7Ti^(0; 0) = which imphes that the photon propagator behaves as free one near 
= 0. Choosing 6Z3 to cancel entire O(e^) correction we find 



SZ^ASS 



lim i7(g2) = 
g2^0 ^ ' 127r2 



— I- 7_E — Zn(47r) + In 
e 



m 



(148) 



and renormalized polarization function in mass-shell renormalization prescrip- 
tion satisfies well known DR: 



00 

nn{0;q^) = n{q^) - lim n{q^) = [ duJ 



(q"^ — u! + ie)ijj 



(149) 



where the absorptive part reads: 



TTp{uj) 



OIQED 



(1 + 2mVt^)\/l - 4mVa;0(u; - Am'^) . 



(150) 



Recall that the one loop nji{0,q'^) represents also self-energy calculated in the 
popular MS scheme for the special choice of t'Hooft scale pt'H = Clearly in 
off-shell momentum subtraction renormalization scheme: SZ^ = iI(C) and the 
DR ((TlUl) is generalized to ((Tin|) . 

In our approximation the fermion propagator entering the photon SDE is 
dressed. Substituting SR for fermion propagator 5* into the photon polarization 

(IHHl) 



[YS{l)j'^Sil-q)] , 



(151) 



and performing the standard steps we can arrive at DR where the weight function 
is now: 



127r2 



da db 



a + b b — a f _ b — a 
1 + ^ ^ 1 + ^^ 



ay{a)a^{b)0 (q^ - {^/^ + Vbf 



(152) 



Clearly form the expression (|152j) follows that the functions P'y,cr^ are nonzero 
only for q^ > 4m^. 



4.10 Numerical comparison 

Like in the ladder approximation considered previously the resulting set of Uni- 
tary Equations have been solved by numerical iterations. The cutoffs yl| = 
Ajj- = 10^M2(0) have been used and the renormalized mass is chosen such that 
M(0) = 1. 
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Figure 9. a) Top figure shows charge renormalization function H = [1 — Ur)'^ for timelike 
momenta for the renormalized couplings: a(0) — 0.1,0.2,0.4. The down oriented peak corre- 
sponds to the threshold 4m^. b) Bottom figure displays charge renormalization functions for 
spacelike momenta as obtained by the solution of the SDEs in Euclidean (EU) and Minkowski 
(S) space are shown in the right figure. Dotted lines stand for one loop perturbation theory. 



The solutions for the running couphngs for timehke momenta as have been 
obtained from the Unitary Equations are displayed in Fig. El For better visual- 
ization the so called photon renormalization function Hn = [1 — 77/j((^; g^)]"^ is 
actually plotted. The spectral and Euclidean solutions are compared in Fig. 
too. Fermion mass function obtained from SDEs is displayed for spacelike regime 
of momenta in Fig. llUI The observed deviation of Euclidean results from the spec- 
tral ones can be explained as a consequence of different cutoffs concepts used in 
the Euclidean and Minkowski treatment. The timelike solution for M is shown 
in Fig. ITIlk ) where the absorptive part is also added. We compare, too, with 
one loop perturbation theory. In this case we use the pole mass obtained from 
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Figure 10. a) Upper panel displays the mass function M{p^) for timelike momenta. Note, the 
sign of the imaginary part, i.e. of (Tvpsip^)). b) The lower panel displays comparison of mass 
function M(p'^) calculated in Euclidean (EU) and in Minkowski (S) formalism for spacelike 
momenta now. 



Unitary Equations as an input for one loop on-shell renormalized perturbation 
theory expression. 

The large growth of the running coupling with its derivative close to the 
cutoff (see Fig. signals a good evidence for the Landau singularity [HH] and 
hence QED triviality. As the consequence of QED triviality the collapse of the 
both Spectral and Euclidean solution is observed. Contrary to quenched approx- 
imation, using the spectral approach here, we have always identified the physical 
pole mass which is probably the consequence of the smallness of the coupling at 
low q^. 
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+ 




Figure 11. Diagrammatical representation of the Schwinger-Dyson equation for gluon self- 
energy. All internal propagators are dressed, as well as the vertices labeled by blobs. Wavy lines 
represent gluons, dashed lines represent ghosts and the solid lines stand for quarks. 

5 Summary on the Unitary Equations 

6 Analyticity of Green's functions in QCD 

6.1 Motivation 

Quantum Chromodynamics (QCD) is the only experimentally studied strongly 
interacting relativistic quantum field theory. This non-Abelian gauge theory with 
a gauge group SU{3) has many interesting properties. The dynamical breaking 
of chiral symmetry explains why the pions are light, identifying them with the 
pseudo-Goldstone bosons associated with the symmetry breaking of the group 
SU{2)l X SU{2)r to SU{2)v (in flavor space). Asymptotic freedom [Ml [Ml 
implies that the coupling constant of the strong interaction decreases in the 
ultraviolet region. For less than 33/2 quark flavors, QCD at high energy be- 
comes predictable by the perturbation theory. However, in the infrared region 
the perturbation theory does not work and nonperturbative techniques have to 
be applied. 

One of the most straightforward nonperturbative approaches is a solution of 
the SDEs for QCD. The extensive studies were undertaken for a quark SDE, 
based on various model assumptions for a gluon propagator. These approximate 
solutions, often accompanied by a solution of the fermion-antifermion Bethe- 
Salpeter equation for meson states, have become an efficient tool for studies of 
many nonperturbative problems, e.g., the chiral symmetry breaking, low energy 
electroweak hadron form factors, strong form factors of exclusive processes, etc. 
(see reviews H El 0^07, 108 and also recent papers [Tl lTMl 111111 ITTTl ITT^ ^ . 

However, to take gluons into account consistently is much more difficult than 
to solve the quark SDE alone. The SDE for the gluon propagator (see Fig. 
is more nonlinear than the quark one. Moreover, the Fadeev-Popov ghosts have 
to be included pJ3J in a class of Lorentz gauges. In a class of covariant linear 
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gauges the contributions to the gluonic polarization function are sketched in 
Fig. 111! Although it is not so difficult to derive and write down the formal 
dispersion relation for vacuum polarization function (especially when one use 
bare vertices) we were not yet able to succeed with a numerical solution of 
Unitary Equations. Therefore, after a brief overview of selected representative 
SDEs solutions in Yang-Mills theory, instead of solving the SDEs we use the 
generalized (nonpositive) spectral representation to fit the spectral function to 
Euclidean solutions obtained in recent lattice simulations. This is a first step 
towards the understanding of analytical structure of infrared QCD. 



6.2 Cornwall's approach 

In this section we mention a symmetry preserving gauge invariant solution to 
gluon propagator obtained by Cornwall two decades ago j^. To our knowledge 
this is the best published example, in which the behavior of the QCD Green 
function in the whole range of Minkowski formalism is addressed within the 
framework of the SDEs and spectral technique. To this end, the gauge invari- 
ant GFs have been constructed by means of the Pinch Technique. To achieve 
this particular spectral Ansatz for GFs were employed. The resulting vacuum 
polarization function is then by construction gauge-fixing independent. We leave 
aside the consistent nonperturbative proof of this statement, but we refer to the 
appropriate perturbative proof in ref. [Ml- Let us now describe the basic steps 
of Cornwall's procedure. 

As in the usual approach, the gauge fixing is imposed in the first step: The 
light-cone axial gauge is employed. Avoiding thus SR for ghosts, only the follow- 
ing SR for the gluon propagator is assumed: 

Gf^AQ) = G{q^) = P^, [ du , ''^'^\ . , (153) 

J - uj + te 

P^iv = -9fjLu + -^, (154) 

[n ■ q) 

where for it holds n^n^ = 0, n • ^ = 0. G{q^) is the gauge-fixing independent 
scalar part of the gauge- invariant propagator G^y{q). The Ward identity 

KPcmO^iM, h) = G^}^{k2) - G^^ih) . (155) 

between the gluon propagator and the three-gluon vertex Pa/syiki, k2, k^) has 
to be satisfied. To enforce this identity the Gauge Technique was used, which 
makes the following Ansatz for the longitudinal part of the untruncated three- 
gluon proper GF: 



G„^7 = Gf"^' {k2)r^f,,^,G^''' {ks 



duj a{uj) 



,2 , ■ Paa' yi^i^ ^2, ks) 2 
K2 — 1^ + «e — uj + le 



(156) 



Recall that the Gauge Technique was originally developed by Salam jll4j 
as an attempt to solve a vector electrodynamics and later extended and ap- 
plied to the scalar IHll and fermion [115^ QED, for later studies based on the 
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Gauge Technique see jll6( \117\ I118( I119j . Substituting the Ansatz (|15f)|) into the 
approximate pinched gluon SDE pT] yields for the inverse of G{q^): 



+ 



q 

ibg 



llvr^ 



1 + 



dw 



[{k + q? 



[{k + qY - 



d^kGik"^ 



where b is the lowest-order coefficient in the (3 function 

33 



(5 = -bg-' + . . . , 
and is the coupling renormalization constant: 

= g^gl ■ 



487r2 ' 



(157) 
(158) 

(159) 
(160) 



Equation ()157() was not solved directly in Minkowski space, it was converted 
into the Euclidean space. Defining 



d{q^) = -G{-q^ 
it was found that the trial propagator 



{161) 



[q^ +m^{q^)]bg^ In 
In 



q^ + Am?{q^ 



m 



In 



4m,2 



12/11 



(162) 



where g, m and A are fitted parameters, is an excellent fit to the numerical 
solution of Eq. H157() . The scale is fixed by the condition d~^{q'^ = C) = C) where 
is the renormalization scale 1(311 . 
Let us stress the main physical results obtained from the combination of the 
spectral Gauge and Pinch techniques. First of all it seems possible to reorganize 
the SDEs of QCD into a form which allows to deal directly with a gauge inde- 
pendent GFs. The second important feature is that the solution has a mass gap. 
There is no singularity corresponding with massless excitation since the gluon 
propagator is finite at zero momenta. For low the mass m{q^) varies within 
500 lb 200 MeV, for large q^ it vanishes, in agreement with the expected (pertur- 
bative) UV behavior. Note that similar behavior was obtained by the operator 
product expansion method several years later in !120j : 
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where Q is the gluon Euchdean momentum, < ^F^^F^'^ > is the gluon conden- 
sate UHB and iVc = 3. 

Let us recall also some recent results supporting the idea (suggested long 
time ago |122| 'l of the dynamical gluon mass. It should be emphasized that the 
gluon "mass" is not a directly measurable quantity, but it can be related to 
other physical quantities. The lattice results P2S1 El dSl 1121 dZl CDI seem 
to support this concept. Parisi and Petronzio |122j were first to introduce it into 
description of J/ip decays. 

6.3 Gluons and ghost in linear covariant gauges 

The solution of SDEs in strong coupling non-Abelian theory like infrared QCD is 
strongly dependent on the gauge fixing of the QCD action and on the truncation 
of equations system as well. Since the coupling increases as one goes from the 
ultraviolet to the lower square of momenta, various approximations can lead to 
the different analytical properties of Green's functions. A quite selfcontaining 
history of SDEs approach and phenomenology in QCD can be find in the paper 
jl28j where the present knowledge about the infrared behaviour of QCD Green's 
functions is summarized in the review (see eg. the Chapter 5). Here we will briefly 
mention the main results obtained in the Landau gauge, noting that the most of 
the recent results are obtained in this gauge. 

In the initial study of SDEs |129j Mandelstam provisionally neglected the ef- 
fect of ghosts. As the consequence of this first approximation it was shown that 
gluon propagator is more singular then simple pole in this case. Using the simi- 
lar truncation this behaviour was confirmed in the studies 113011131] . It was first 
suggested in |132j that such behaviour is not an inherent property of nonpertu- 
bative dynamics of Landau gauged QCD but rather of deficiency of truncation of 
SDEs system. It was suggested that this is not gluon, but rather ghost propaga- 
tor that is highly singular in the zero momentum limit. The detailed analytical 
and numerical study of the coupled ghost and gluon equations in Landau gauge 
with bare vertices was made (133 , 1341. It was actually confirmed that the ghost 
loop dominates the infrared behaviour of the gluon propagator. In recent papers 
[3 ini 11361 11371 1138j studies of the coupled SDEs for gluon and ghost propagators 
in the usual Landau gauge in Euclidean space were performed in various ap- 
proximations. An interesting result of the investigations in Landau gauge is the 
observation, that there is no qualitative difference of the solutions found with 
bare vertices or with vertices dressed by the use of Slavnov- Taylor identities. Note 
that the similar behaviour was observed in the wide class of ghost-antighost sym- 
metric gauges jl39j . The attempt to estimate the effect of two loop irreducible 
contributions was made in [^. The numerical result indicates slight changes in 
the intermediate regime. Note that the inclination of this result towards the 
lattice data was shown already in the introduction of presented paper. 

6.4 Analytical continuation of Euclidean solutions 

The assumption of GF's analyticity and hence the assumption of the propa- 
gator spectral representation was key point for the formulation of the Unitary 
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Equations and the solution of SDEs in Minkowski space. The requirement of the 
propagator (or invariant charge) analyticity is supposed be a good assumption 
in QCD for many reasons jl4()[ I14l1 1142| \T^^ and the analytical structure of the 
GFs obtained form Euclidean SDEs in QCD is under current investigation j76j . 

In this Section we perform analytical continuation of recent unquenched lat- 
tice results [Ij for gluon propagator to the Minkowski space. Note that the similar 
fit has been already performed for the case of pure gluodynamics [33 where the 
quenched lattice date ^Hl has been used. In the paper P the authors considered 
three flavor QCD. For this purpose one massless and one massive quark have been 
taken into account. The main difference that follows from unquenching of lattice 
gluon propagator is the partial suppression of infrared enhancement plateau. For 
some details see the original paper where the comparison is made. Of course, 
taking the effects of quark-antiquark pairs into account it has no influence on 
the main characteristic feature of the Landau gauge gluon propagator: it is ob- 
viously infrared vanishing function. Therefore, in contrast to the propagator of 
a stable unconfined particle, we assume that the spectral function is a smooth 
real function, i.e., propagator has not simple pole structure. 

In fact, what we actually do, it is the continuation of the data obtained for 
a real negative momentum semiaxis to the border of analyticity domain. This is 
the standard situation when one make prediction from experimentally obtained 
data to the experimentally inaccessible regime (the extraction of nucleon spectral 
function from the electromagnetic form factor is a good example) . This approach 
has been already extended to the case of gluon spectral function of SU (2) Yang- 
Mills theory ref. I I."! . However, in our case, the smallness of errors of the lattice 
simulations and the value of lattice spacing as well allows to use the data directly 
for our purpose. No extrapolation to higher square of momenta is made and we 
neglect the computational errors as well. 

The gluon spectral function cr(uj) and the most probable analytical fit is 
obtained from the minimization of 



Ifit 



^[H{xi)-H{xi)] , (164) 

i 

(165) 



where H is the gluon form factor defined as H{x) = xG{x) , H{x) are the 
appropriate lattice data evaluated for the square of momenta x = while H{x) 
is obtained through the fit of the spectral function which enters the generic SR 
for the gluon propagator in the Landau gauge which reads 



Gab^q) = Sab 



-9-' + 'X 



G{q') = I dc. , ^^"^^ . , (166) 
' q'' — LO + ie 







where A, B are color indices. The Ansatz ()166() should be considered as the 
generalized spectral representation, since we do not assume (and do not obtain) 
the spectral function (t{uj) positive for all values of a;. For completeness let us note 
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Figure 12. The extrapolated gluon form factor H(x) — xG{x) at timelike momentum axis. 



that we use unrenormalized data, i.e. the suitable but for us irrelevant constant 
prefactor should be added when considered for some application. 

To achieve a reasonable accuracy we use the following iterative minimizations 
procedure: The function Hj evaluated at some J-step of our iteration cycle is in 
each case obtained by the previously obtained fit for spectral function ,i-e. (Jj-i 
such hat the relation for Hj reads: 



Hj{xi 



Xi 



duj 



Xi+UJ 



1=0 



u> + c 



(167) 



where are dimensionless constant coefficients of Legendre polynomial which are 
fitted during the minimization procedure, c is positive massive parameter which is 
use for the fit as well. Some suitable zeroth iteration input is used for ctq. For this 
purpose we use the convenient choice of the form: ao{iv) = K ^°^ ui+A^^^ ^ with 
some suitable parameters 7, A. Assuming that the procedure is convergent 
then the best fit is then obtained when oq = 1 and other coefficients of Legendre 
polynomial are zero, which should be true irrespective of the value of the constant 
c. We actually observe the relatively fast convergence of the described process 
and we stopped our minimization procedure when oq = 1.0006 and other a's 
are at most few of 10~^. Such minimization procedures is much powerful then 
to perform minimization at once time with some 'more suitable' function even 
with larger number of the fitting parameters. As a bonus the computer time is 
reasonable short (seconds for office PC). 

For a vanishing gluon propagator we have the following "sum rule" for a 



G{0) 



a{uj) 



UJ 



~ 0. 



(168) 
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Figure 13. Gluon propagator for Nf = 3 QCD, the lattice data pOj are represented by dia- 
monds. The fit (at spacelike) and the extrapolation (at timelike) is labeled by solid line. 



which is clearly impossible for the positive a. Note that the violation of positivity 
was confirmed by some other lattice simulations jl44j , where it was indicated by 
the study of Schwinger function. The extrapolated solution is plotted in Fig. [T2l 
The spectral function has a positive peak at the infrared and becomes negative for 
asymptotically large p^. The infrared behaviour describes the shortliving gluon 
excitation and is consistent with the absence of gluon in the physical spectrum 
while the large momentum behaviour is the feature already expected from the 
perturbation theory. Clearly, the observed shape of the gluon propagator in Lan- 
dau gauge is in accord with our physical expectations, e.g. with the confinement 
of color object and asymptotic freedom as well. 

In Fig. El we also display our fit at the spacelike regime and provide the 
comparison with the original lattice data ^1 . The real part of extrapolated gluon 
propagator is added for comparison too (negative semiaxis now). From this 
we can see that we were not able to fit the lattice data exactly. The origin of this 
small discrepancy can be a combination of several effects: it can be artifact of 
our iteration/minimization procedure, it can signal the the complex singularities 
however positioned rather close the real positive momentum square semiaxis, 
partially it can be a consequence of finiteness of the lattice spacing and finite 
lattice volume. This complications could be clarified in future studies but they 
have a little significance for the main statement: The Landau gauge gluon prop- 
agator is the analytical function in the whole complex plane with the exception 
of (or small vicinity of) a real positive semiaxis of square of momenta and it can 
be parameterized (or at least reasonable approximated) by nonpositive spectral 
representation. 
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7 Summary, conclusion and outlook 

An important question, in field theory in general, and in the SDE approach 
particular, is the connection between Euclidean-space and physical Minkowski 
space. As we explicitly demonstrated it is clear that to perform SDEs studies 
directly in Minkowski space is possible. 

We began with an introduction to the analytical concept in the formalism 
of Schwinger-Dyson equations with specific emphasizes to the derivation and so- 
lution of the Unitary Equations in scalar models and QED. For the later case 
we presented studies of strong-coupling QED in different approximations. One 
of them was unquenched QED, which represents the study where the running of 
the coupling has been correctly taken into account. In this case the QED trivi- 
ality plays its own crucial role in the asymptoticly large spacelike momenta and 
we explain how to deal with this trivial theory correctly within the formalism 
of spectral representations and dispersion relations. The Unitary Equations pre- 
sented are systematically improvable by considering further and further skeleton 
contributions to the effective action. If the analyticity is correctly maintained 
then the all momentum integrations can be performed analytically. For instance, 
let us consider the two-loop improved approximation of our scalar model con- 
sidered in the Section 3. In this case it is sufficient to add the contribution from 
the two-loop graphs to the weight function (|58|) . i.e. 

Pip")^ p{p')+Psun{p') + ..., (169) 

where Psunip^) is given by rhs. of Eq. pH]) and the dots stand for the terms 
stemming from the two other known skeletons not explicitly considered here. 
This is a clear technical advantage of the presented method: we can immediately 
apply the known original results for absorptive parts of Feynman diagrams when 
they were already evaluated in perturbation theory. 

For an explicitly massive theory and up to the certain values of the coupling 
constant the accurate comparison of SDEs solutions at spacelike momentum 
axis with the solutions performed by more conventional strategy in Euclidean 
space gives us strong belief in the relevance of presented methods. In this cases 
the only limitation of spectral method is given by the increasing number of 
spectral integrals in the expression for DR that is directly dictated by the number 
of internal propagators. Nevertheless, within contemporary power of computer 
facilities many of exciting calculations can be made and remain to be done. 

The strong coupling QED is often regarded as an reasonable pedagogical tool 
for SDE studies and their application to QCD. We also attempted to extend the 
discussion to QCD where the direct Minkowski space formulation and solution 
is presently prevented by the partial lack of information on the timelike axis. 
While many important steps have been accomplished, it is apparent that much 
more needs to be done. Furthermore, at present time and up to a few exeptions, 
the vast amount of studies dealing with the deep nonperturbative phenomena 
like chiral symmetry breaking and dynamical mas generation were carried out 
in Euclidean space. Regarding the real QCD, it is not clear when or if it will be 
possible to perform such symmetry breaking SDEs studies directly in Minkowski 
space. 
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Another interesting prospective is to further explore the timehke infrared be- 
havior of Green's functions by analytical continuation of presented lattice studies. 
We offered the solution for Landau gauge gluon propagator, but the information 
about the timelike structure of quark propagator is currently unavailable. The 
present day great improvements in the QCD lattice studies could provide signif- 
icant benefits at this area in the recent and in the near future. 

We also collected the necessary information and illustrated how it can be used 
in relativistic bound state calculations using the spectral representation for solv- 
ing the Bethe-Salpeter equation in (3+1) Minkowski space. To further develop 
the method, it would be interesting to extend it to more complicated kernel: 
trying to include the cross boxed contributions, s and u channel interactions etc. 
It is already known that the "spectral approach" used here is suitable even for 
more complicated systems. The same or very similar tricks and technology when 
successfully performed in the Minkowski calculation could become great goals in 
hadronic physics in the not too distant future. As it is already known the obsta- 
cles due to fermionic degrees of freedom can be overcome. The solution of the 
spinless bound state (pion) composed from two quarks interacting via dressed 
gluon exchange is under auspicious consideration. 
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Appendix A: 

The standard trick is to re- write the integral J(p^) in terms of the Feynman parameterization: 



(2nY J {2iiY Jo [kl + klx{l-x)~a^x-a2{l-x)+ief 
rf'*fc2 /■'^^ iD{k2+p,a3)r{2-d/2) 



{'^■^Y Jo (47r)d/2[x(l - x)]2~d/2 ^fc2 _ J2J_ _ £^ + jgj ' ""^^ 

d'^ki i{4-K)-''/^r{3-d/2)y^-''^^ 



dxdy 



{2-Ky - j:)]2-d/2 

1 



[fc2 + p2j^(i _ J/) _ [ Jii. _ - 03(1 - y) + le 

(_)2-d/2j.(3 _ ^)yl-d/2 



3-d/2 



dxdy 



(47r)''[s(l - 3;)]2-'«/2[y(l - y)]3-d (p2 _ ^ + jg)^ 



^^Sixl^a, . s{x) = ^ + ^. (A.l) 
1 — y y 1 — X X 



To make the subtractions (I28II it is convenient to use the following Feynman formula 

qjB - A) 
' [{A - B)w + B]i+'^ 



-B--=l' . (A.2) 
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Using variable z for the mass subtraction and the variable u for the subtraction of the term 
with derivative, we can immediately write down a formula for renormalized J: 



J«(P^) = J(P^)-J(C)-^ 



(p' - C) 



dxdydz dur{5 - d) z{-\f-''l'^ y'-'^l'' {p" - 0' 



(47r)<* [x(l - 2;)]2-d/2 [y{l - y)]3-d [(p2 _ Q^y^ + ( - ^2 + is]' 



(A.3) 



In four dimensional spacetime the integral l|A.3^ is already finite and the limit d ^ 4 can be 
simply taken. 

In addition we make the following substitution: lj — + ^^-^ . After changing ordering of 
the integrations du <^ dz and integrating over z gives: 



Jr{p ) = / dxdy 



dio 



(47r)4(a;-C)2(p2-t^ + ie) 



(A.4) 



Notice that all x, y dependence is now hidden in the lower bound Q, which is defined in last 
line of IIA.lll . It is still possible to evaluate one more integral analytically, we have chosen to 
take the integral over y. Changing again the order of integrations we get for boundaries: 



dy 



dio - 



duo 



where j/+ are the roots of the equation ui — Q = G: 



dy 



y± = 



uj — s + a-i ± ^ X{S {x) , uj , a-i) 



2ij 



(A.5) 



(A.6) 



where X(x, y, z) is the triangle function defined in 12411 . After the explicit integration over y and 
changing order of the last two integrals we arrive at the following result: 



du ■ 



ip^-Cfpj{oj) 



1 



(A.7) 



(47r)4 
as In 



dx 



{cj - S{x) +a3)^\iS{x),Lo,a3) 



tjj — S{x) + Qfs + \/\(S{x),uj,a'i) 
LO — S{x) + a:i — ^ \{S{x),LU,a3) 



(A.8) 



x± 



If the omitted prefactor is restored Eq. IIA.711 we can simply identify the one dimensional integral 
representation for imaginary part of IIr/tt. 

This result can be rewritten into the form already available in literature |24 | .|25 | . To this 
end we have to make the following substitution x ^ s such that 

s = S{x) ■ ^ 



s + ai-a2 + ^A(a2,ai, s) — — 2x 
for X > x{s = ( V"i + VQ?2) ) 



2s 



^ ^ 2s 
which leads to the relation given in |24 | .|25 | : 

pj(i^) = . . . , / ds 



for X < x{s = (a/oi + 1/02) 



(47r)4 



^X{s,LO, Q;3)^A(q2, S, Qi) 



© U - (E 



1/2x2 



(A.9) 



(A.IO) 
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